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DYNAMIC HEAVE-PITCH ANALYSIS OF 
AIR CUSHION LANDING SYSTEMS 

By K. M. Captain, A. B. Boghani, and D. N. Wormley 
Foster-Miller Associates, Inc. 


1. Introduction 


As part of the effort to advance Air Cushion Landing System 
(ACLS) technology, NASA has initiated a program to develop analytical 
tools to help evaluate ACLS dynamic performance. This report 
describes the first two phases of this program, which are now complete. 

The objective of these phases was to formulate a fundamental 
analysis of the dynamic behavior of the ACLS and develop a computer 
program to carry out the dynamic simulation. First, the heave (vertical) 
motion of the ACLS was analyzed, and the analysis was then extended, 
and a coupled heave-pitch model was formulated. 

The mathematical models are based on a fundamental analysis of 
the body dynamics and fluid mechanics of the aircraft- cushion- runway 
interaction. The analysis takes into account the air source character- 
istics (fan, etc. ), flow losses in the feeding ducts, trunk and cushion, 
the effects of fluid compressibility, and dynamic trunk deflections, 
including ground contact. The computer program developed is capable 
of simulating the dynamic motion of an ACLS-equipped aircraft caused 
by landing impact and taxi over an irregular runway, using input data 
such as cushion and trunk geometry, aircraft weight, fan characteristics, 
runway surface profile, etc. 

The program can be used in three principal ways: 

1. To determine static ACLS characteristics (equilibrium 

height, stiffness, static pressures, etc. ), aid in fan 


selection, and determine allowable limits for equilibrium 
cushion loading. 

2. To evaluate dynamic landing and taxiing performance 
(g loading, heave and pitch motion, trunk deflections, 
hard surface clearance, etc. ), including the vibration 
caused by runway irregularities. 

3. To determine optimum values of design parameters 

(e. g. , hole sizes and configuration, trunk shape, etc. ) 
for improved dynamic performance (i. e. , design guide- 
lines ). 

The types of performance results that can be obtained from the 
computer program, which is described in the Appendices, are shown in 
Table I. Illustrative simulations of a scale model ACLS have been 
carried out and are presented in Section 3. These simulations show the 
results obtained from the program for three typical cases of interest 
a zero speed heave drop test, a 22. 4 m/s (50 mph) landing impact, and 
taxi over an irregular runway. 

During the next phase of this program, the analysis will be 
extended to include coupled heave-pitch- roll simulations, and the anal- 
ytical models will be verified and refined based on results obtained with 
a test cushion at NASA- Langley. After model verification, a series of 
additional simulations are planned to investigate a variety of potentially 
attractive ACLS configurations and develop guidelines for improved 
designs. 
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Simulation Capabilities 
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2 . 


Analysis 


The analysis outlined herein is the generalized pitch-heave 
analysis, further details of which are given in Appendix C. The 
pure heave analysis - a special case of the heave-pitch analysis - 
obtained by setting the torque and angular motion terms to zero- 


2. 1 Basic Configuration 

The basic ACTS configuration analyzed is shown in 
Figure 1. The model includes four primary subsystems - i. e. , the 
fan, the feeding system, the trunk and the cushion. The configuration 
of these systems has been chosen sufficiently general so that they can 
represent a wide variety of practical designs. Air from the fan flows 
through the ducts and plenum (feeding system) and enters the trunk. 
The trunk has several rows of orifices that communicate with the 
cushion and atmosphere. Thus, the airflow from the trunk has two 
components - one part entering the cushion and the other leaking 
directly to the atmosphere. The cushion flow exhausts to the atmos- 
phere through the clearance gap formed between the trunk and ground. 
In addition to the basic flows described above, two other flows have 
been included in the model, for generality. These are the plenum 
bleed flow and the direct cushion flow. Plenum bleeding causes some 
of the air to flow directly from plenum (fan outlet) to atmosphere, 
and has been used in some designs to improve the dynamic character- 
istics of the air supply system. Direct flow from the plenum to the 
cushion can also improve dynamic response. 


In plan, the cushion has an oval shape, made up of a 
rectangular section with semicircular ends. a and b are the horizontal 
and vertical distances between the points of attachment of the trunk to 
the aircraft body. The initial (undeformed) trunk shape is defined in 
terms of the above two parameters, and the perimeter f and height h 
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(a) Plan View 

Figure 1. Basic ACLS Configuration 


Airflow from Fan 



(b) Front View in Cross-Section (Section A-A in Figure 1(a)) 
Figure 1 (Concluded), Basic ACLS Configuration 





as shown. is the (uniform) spacing between the rows of peripher- 

ally distributed orifices. The number and orientation of the orifices 
can be selected independently in terms of the number of orifice rows 
(N^), the number of orifices per row* (N^), and the orientation para- 
meter (f ). The cushion volume consists of two parts: an active 

(dynamically varying) region and a dead (static) region. The active 
cushion^^olume depends only on the trunk shape and ground profile, and 
is computed by the program from cushion geometry. The dead volume 
(shown in Figure 1) includes recesses in the cushion cavity, and is a 
design variable. 

2. 2 Assumptions 

Order- of- magnitude analyses and available test data have 
provided the initial basis for determining the principal assumptions of 
the heave-pitch analysis. These assumptions are summarized below. 

(a) Computation of the trunk and cushion parameters 

(trunk volume, cushion area, etc. ) is carried 
out by dividing the trunk and cushion into seg- 
ments, as shown in Figure 2. The parameters 
are first calculated for each individual segment, 
and then added together to obtain parameter 
values for the full cushion. 

The ground under any particular segment is 
considered parallel to the hard surface, and at 
an elevation corresponding to the ground profile 
at the segment center projection on the reference 
plane (as shown in Figure 3). This assumption 
represents the ground surface and the hard 
structure of the cushion by a series of short, 
parallel sections which, when chosen sufficiently 


7 









small, closely approximate the actual ground 
profile and hard surface orientation. 

The two types of segments are shown in Figure 2: 
Rectangular segments in the straight portion of 
the cushion, and pie- shaped sections at the 
curved ends. 

The trunk and cushion pressure force components 
for each segment are found from the products of 
the appropriate pressures and areas, and are 
represented by concentrated forces acting at the 
respective centers of pressure, as shown in 
Figure 4. 

(b) The flow analysis is based on a lumped para- 
meter model of the ACLS as shown in Figure 5. 
Plenum, trunk and cushion pressures are assumed 
uniform (though unequal and dynamically varying). 
The plenum, trunk and cushion cavities are 
represented by their capacitance (volume). Pres- 
sure losses in the ducts and in the entrance and 
exit regions of the chambers are represented in 
terms of lumped orifice resistances. 

(c) For typical ACLS designs, the fractional pressure 
drop across the orifices (Ap/p) is small (usually 
less than 0. 1), and changes in air density across 
the orifices will be negligible. Therefore, the 
orifice flow Q can be found from the incompres- 
sible flow quadratic relationship 




: Sin0 2 

runk Contact with Ground 






where 



C d 

- discharge coefficient 


A 

- orifice area 


P 

- mean air density 

and 

Ap 

- pressure drop across orifice 


Within the chambers, however, air compress- 
ibility cannot be neglected, because dynamic 
density changes (dp/dt) will be significant. 
Therefore the effects of density changes in 
the plenum, trunk and cushion are included 
in the analysis. Density changes are deter- 
mined from pressure changes through the 
polytropic relationship p/p = constant, where 
the exponent K lies between 1 (isothermal 
expansion) and 1. 4 (adiabatic expansion). 

(d) The trunk is modeled as a massless unstretch- 

able membrane capable of bending freely. Initial 
calculations carried out for selected trunks of 
current interest indicate that the deforming 
pressure forces are very large compared to 
the inertia of the trunk material, so that 
changes in trunk shape will occur almost 
instantaneously, and the massless approxima- 
tion will be valid. 

With the assumption of no trunk stretch, the 
trunk length around the cushion periphery will 
be constant. This means that, for uniform 
motion, every trunk element will remain in the 
same lateral position, since any lateral trunk 

motion would require peripheral stretching of 

*For elastic trunks, the inelastic "frozen" trunk model described, requires 
modifications to include elastic effects. 
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the trunk membrane, Therefore, to a first 
approximation, the shape of the trunk cross 
section, when out of ground contact, is 
'’frozen 11 (i. e. , independent of the pressure) 
and depends only on the initial prefabricated 
trunk shape (Figure 6a). When ground contact 
occurs (Figure 6b), the trunk material in the 
contact zone conforms with the ground surface 
by crumpling, while the part of the trunk not 
touching the ground remains undeformed. 

Initial observations of the deformation charac- 
teristics of the two trunks cited earlier support 
the idealized model of trunk behavior described 
above. 

(e) The air source is characterized in terms of a 

static pressure-flow relationship, with hysteresis 
losses included to represent the effects of stall. 
The general source characteristic is shown in 
Figure 7. Curve AB represents the normal 
(unstalled) operating regime, while curve CD 
represents the stalled characteristic. The shape 
of the curves depends on the type of the air 
source, and by selecting an appropriate stall 
point A, recovery point C and curve shapes AB 
and CD, a variety of stalling and non- stalling* 
air sources including axial and centrifugal fans 
can be simulated. 

For any pressure, the flow is found by using 
the appropriate (unstalled or stalled) character- 
istic. At the start of the simulation (stall-free 


When point C coincides with point A, stall is suppressed. 
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General Air Source Characteristic • 





initial conditions), the appropriate pressure-flow 
relationship is given by curve AB. When the 
pressure exceeds the stall pressure QP1, the flow 
decreases suddenly, and it is found from the 
stalled characteristic CD. Stalled operation along 
characteristic CD continues as long as the pres- 
sure is above the recovery pressure QPZ. When 
the pressure drops below QPZ, the flow increases 
suddenly (i. e. , recovery), and the pressure-flow 
relationship is given by the unstalled characteris- 
tic AB. The above discussion indicates that down- 
stream pressure variations large enough to cause 
stall and recovery result in a net energy loss due 
to hysteresis (see Figure 7). 

The present analysis is based on the initial 
assumption that the fan flow changes simultaneously 
with pressure. In practice, however, the effects 
of fluid inertance will introduce lags in the flow, 
particularly during the stall and recovery transi- 
tions. The effects of these lags will be to slow- 
down fan stall and recovery, and hence slow down 
passage around the hysteresis loop shown in Fig- 
ure 7. However, since typical fan flow lags are 
estimated to be small compared to the character- 
istic periods of ACLS motion, the lags will have 
only a small affect on the predictions of overall 
landing dynamics and aircraft g loading. Subse- 
quently, detailed stall investigations may require 
a more advanced model which includes fluid inert- 
ance and flow lags. The basis for development 
of this improved model will be established through 
dynamic fan tests scheduled later in the program. 
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<f) 


Five mechanisms of energy dissipation are 
included in the analysis. 


(i) Fan stall and recovery losses (see 
above). 

jjj 

(ii) Aerodynamic drag of the cushion. A 
square law relationship is assumed, such that 
the drag force F is given by 


F = C D A plP v2 


where 


C 


D 


A 

P 


P 

and 

V 


heave drag coefficient 

projected area on which Cjj is 
defined 

ambient air density 
heave velocity cushion. 


(iii) Damping due to trunk crumpling during 
ground contact (see Figure 6). In this case, 
the damping force F is assumed to be linearly 
proportional to the trunk segment deformation 
velocity, V . The trunk damping force acting 
on the aircraft is thus given by 



Drag relationships are preliminary and primarily valid for zero speed 
drops. In subsequent phases, more detailed models of the aerodynamic 
characteristics are planned for inclusion in the model. 
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where B is the damping constant for each trunk 

z i 

segment, and the summation is jcarried, out over 

all the segments. The damping' constant is 

estimated for the trunk sizes and configurations 

of interest by dimensional analysis, using test 

data obtained with prototype cushions. The 

trunk damping force also develops a torque 

around the CG. 

(iv) Energy losses in the orifices. 

(v) Friction losses due to trunk-ground 
contact. The friction force which arises at 
the trunk-ground interface results in a hori- 
zontal retarding force and torque at the CG. 

(g) Because of the presence of brake tread material, 

trunk imperfections, ground irregularities, etc. , 
sealing of the trunk orifices and the cushion-to- 
atmosphere exit area will not be complete, even 
when the trunk is in ground contact. The effects 
of incomplete orifice closure are taken into 
account in the analytical model through blockage 
factors that allow some leakage flow to occur 
even when the orifices are nominally closed. 

2. 3 Analytical Development 

The analysis provides 

(a) The relationships that determine the static 

cushion characteristics (pressures, flows, etc. ) 
existing at equilibrium. These relationships 
are also used to determine the initial conditions 
for the simulation. 
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(b) 


2 . 3 . 1 


The differential equations of flow and motion 
(state equations ) from which the pressures, 
flows, displacements, accelerations, etc. can 
be determined as functions of time. 

Static Model 


The equilibrium conditions are found as follows: 

(a) By applying the steady- state flow con- 
tinuity equations to the plenum, trunk 
and cushion cavities (see Figure 5). 


^fan 

Q plat + Q pltk + 

Q pltk = 

^tkch + ^tkat 

^chat 

^plch + ^tkch 


(b) By satisfying the fan flow constraints. 


Q fan = f < P fan> 


i. e. , where the fan flow Q ^ n and pres 

sure rise P, are determined from th< 
fan 

characteristic fan curve. 

(c) From the static force balance equation 
Forcn = <P ch A ch + P^) co.* 
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where 


Forcn 


P 

A 

P 

A 


ch 

ch 

tk 

tkcn 


4 > 


aircraft weight (in equil.) 

cushion pressure 

cushion area 

trunk pressure 

trunk area in ground 
contact 

pitch angle. 


(d) From the static torque balance equation. 


Torn = o= 


2(M+N) 

I 


i=l 


2P ch ( A ch (i l)( X eh (i) - CC ) 


+ 2P tk ( A tke„ (i >) K^- CC ) 


where 


Torn 


A ch (i) - 


A tkcn (i| 


X.hli) - 


**<*> - 


torque about CG (zero 
in equilibrium) 
cushion area correspond- 
ing to ith segment 
- trunk contact area 
corresponding to ith 
segment 

distance between the 
center of pressure of 
the ith segment of the 
cushion and the geometric 
center of the cushion, 
distance between the 
center of pressure of 
the ith segment of the 
trunk and the geometric 
center of the cushion. 
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2 . 3 . 2 


State Equations 


The state equations are derived from the dynamic 
ACLS model (Figure 8) as follows. 

1. Plenum Flow Continuity 

The net inflow equals the rate of increase 
of fluid mass within the plenum 



where p is the mean air density 

2. Trunk Flow Continuity 

Similar to (1) above 


dt 


* PV tP ( Q pltk ' Q tkch " Q tkat) p 


3. Cushion Flow Continuity 

dF (PV ch ) "( Q plch + Q tkch " Q chat) P 

4. Force Balance about the eg 

M a^r Y cg = < P ch A ch + P .k A tkcn> 


- - 7 C D> P ( TI*) 


Aerodynamic Drag 
Component 


Foret 
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Trunk Damping Component 


Pressure, Flow 



Characteristics Damping Force 

and Torque, Ground 
Contact Friction 
Force and Torque 


Figure 8. Dynamic ACLS Model 
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where 


2(M+N) 

Foret = 2 ^ 1 Forct(i) 

i=l 


and 


Forct(i) 


B 


4(M+N) 


(x_<i>-cc)) 


if the ith segment is in 
ground contact 


if the ith segment is 
not in ground contact. 


Torque Balance around the eg 

P c^cM X c^- CC ) 

P tk\kc,p( x «k< l >- Cc ) 
Ground friction torque 

Torque due to Aerodynamic 
Drag force 
“ T orqt 


2 2(M+N) 

Inert • —jr = 2 Y! 

dt i= 1 


2(M+N) 

♦ z E 

i=l 

+ Torf 


Trunk Damping Torque 


where 


Torf = _ 2 


Torqt = 2 

where 


Torqt(i) = 


2(M+N) 

I I VWn (i > •* ( Y gh (i > +GG ) 

2(M+N) 

^ Torqt(i) 

i=l 


4(KRN) (tP + If ( x c*< i >- cc )) 


X (i) CC^ ^ the ith segment 
. tk' 1 ' - » is in ground contact 


0 if the ith segment is not 
in contact 
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3. 


Illustrative Simulations 


A computer program incorporating the heave and heave-pitch 
analysis has been developed. With this program, the dynamic behavior 
of an ACLS-equipped aircraft (g loading, trunk deflection, cushion pres- 
sure, etc. ) can be determined for landing impact and taxi over an 
irregular runway, using input data such as cushion and trunk geometry, 
aircraft weight, fan characteristics, runway surface profile, etc. The 
organization and use of the computer simulation program is described in 
Appendix A. 

Three types of illustrative simulations have been carried out, 
to demonstrate the capabilities of the program. They are 

(a) A drop test simulation, (zero forward speed, pure heave. 
Torque and angular motion terms = 0. ) 

(b) A landing impact simulation. (With forward speed and 
initial angle of attack. ) 

(c) A simulation of aircraft dynamics when crossing a 
runway obstacle. 

In the above simulations, the input parameters corresponded to 
a model cushion that will be tested in a subsequent phase of this pro- 
gram to verify and refine the analytical model. The general geometry 
of the model cushion is defined by Figure 1 and the detailed geometric 
input parameters are listed in Appendix F. The computations have been 
made in English units and converted to SI units. 

3. 1 Drop Test Simulation 

The drop test simulation of the cushion has been carried 
out for a static load of 1220 newtons (275 lbs. ) and a drop height of 0. 152m 
(6 in. ). The corresponding impact velocity is about 1. 5m/sec. (5 ft/sec. ). 
The simulation results are shown in Figures 9 through 13. 
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The static characteristics show that the cushion pressure 
increases with load, and the flow and hard surface clearance decrease 
with load, as expected. The maximum load capacity of the cushion 
(i. e. , the peak load for which stall-free fan operation is possible) is 
about 4000 newtons (900 lbs. ), which is about three times the static 
load. The time history of cushion motion shows that the peak trunk 
deflection is about 38 mm (1. 5 in), which is well within the static hard 
surface clearance of 185 mm (0.611 ft). The period of one cycle of 
oscillation is about 0. 15-0. 2 sec, which corresponds to a characteristic 
heave frequency of about 5-6 hz. The peak acceleration is about 
50 m/s^ (5 g). At impact, the cushion pressure increases to about four 
times its equilibrium value, and this causes the fan to stall. As the 
pressure drops, the fan recovers, and remains in the stall-free operat- 
ing regime throughout the remainder of the simulation. Prolonged heave 
motion excited by repeated fan stall and recovery is thus inhibited. 
Although the impact disturbance begins to die out after the initial bounce 
(i. e. , the system is dynamically stable), the low cushion damping indi- 
cates that several additional cycles will be required before the cushion 
reache s equilibrium. 

3. 2 Landing Impact Simulation 

The landing impact simulation has been carried out 
for a static load of 1220 newtons (275 lbs), and an initial eg 
height of 0.52 m (1.7 ft) (touchdown sink speed of 1.52 m/s). The 
touchdown (forward) speed was chosen at 22.4 m/s (50 mph), with an 
initial angle of attack of 5°. The simulation results are shown in Fig- 
ures 14 through 18. 

The static characteristics (Figure 14) illustrate that the 
eg elevation Increases as the load reduces. The slope of the load- 
deflection curve (stiffness) is smaller for a non-zero pitch angle than 
for a zero pitch angle, because non-uniform trunk contact results in a 
lower restoring force than uniform trunk contact. 
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Static cushion load 1220 newtons (275 lbs) 
0. 152 m (6 inch) level drop 








Acceleration 


Static cushion load 1220 newtons (275 lbs) 
0. 152 m (6 inch) level drop 



Time History of Acceleration 
Figure 11 
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Figure 14. ACLS Static Characteristics 



Static cushion load = 1220 newtons (275 lbs) 
Landing velocity = 22, 4 m/s (50 mph) 
Touchdown sink rate = 1.52 m/s (5 ft/sec) 
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Static cushion load = 1220 newtons (275 lbs) 
Landing velocity = 22.4 m/s (50 mph) 
Touchdown sink rate = 1. 52 m/s (5 ft/sec) 



Figure 16. Acceleration during Landing Impact 




Static cushion load = 1220 newtons (275 lbs) 
Landing velocity = 22.4 m/s (50 mph) 
Touchdown sink rate = 1.52 m/s (5 ft/sec) 
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Figure 17. Cushion Pressure during Landing Impact 



Figure 18 shows the time history of heave-pitch motion 
caused by landing impact. Initially, the cushion has a positive angle 
of attack. As it touches the ground, a clockwise torque acts 
upon it and causes the nose to pitch down. This torque arises because 
the rear of the cushion touches the runway before the front. After 
touchdown, the cushion begins to recover, and the heave and pitch 
motions begin to damp out. 

Figures 16, 17 and 18 show the acceleration, cushion 

pressure and flow during landing. The pressure and acceleration 
build up as the aircraft descends. The increasing cushion pressure, 
however, causes the fan to stall, which reduces its output, and sub- 
sequently decreases the pressure and acceleration. As the pressure 
drops below the stall pressure, the fan recovers and the pressure and 
acceleration build up to a second peak, and then approach their respec- 
tive equilibrium values. 

3. 3 Obstacle Crossing Simulation 

The obstacle crossing simulation has been carried 
out for a static load of 1220 newtons (275 lbs). Prior to 
obstacle impact, the aircraft is assumed to be moving straight and 
level, with a velocity of 22.4 m/s (50 mph). The obstacle is repre- 
sented by a rectangular cleat 0.4 m (1.3 ft) long and 89 mm (3.5 in) 
high. The simulation results are shown in Figures 19 through 22. 

The time history of heave-pitch motion (Figure 19) shows 
that the aircraft (cushion) begins to pitch forward (clockwise) as the 
cushion first impacts the obstacle. This is because the friction force 
due to obstacle contact and the unbalanced (vertical) pressure force 
acting on the rear trunk give rise to a clockwise torque about the eg. 
The entry of the obstacle into the cushion also causes the cushion and 
trunk pressure force components to increase, which results in an 
upward heave motion of the aircraft. The upward motion continues 
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Static cushion load = 1220 newtons (275 lbi 
Taxiing velocity = 22. 4 m/s (50 mph) 



Figure 19. Cushion Motion during Taxi over Irregularit 







Static cushion load = 1220 newtons (275 lbs) 
Taxiing velocity = 22. 4 m/s (50 mph) 
Irregularity = 0. 4 m (1. 3 ft) long, 89 mm 
(3. 5 inch) high cleat 


C 
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until after the trunk leaves the ground. The upward force 
components then vanish, and the aircraft begins to descend. The 
initial pitching torque causes the pitch angle to build up to a maximum 
of about 7°, when the leading edge of the cushion contacts the ground 
and provides the restoring torque that causes the pitch disturbance to 
die out. The heave disturbance also begins to damp out. 

Figures 20, 21 and 22 show the accelerations, cushion pressure 

and flow wl*ile crossing the obstacle. Initially, as the cushion impacts 
the obstacle, the pressure and heave acceleration build up. The 
increasing cushion pressure, however, causes the fan to stall, which 
reduces its output and subsequently decreases the pressure and accelera 
tion. As the pressure drops below the stall pressure, the fan recovers 
and the pressures and accelerations reach another peak at the second 
bounce, and then approach their equilibrium values. 
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4 . 


Conclusion 


The effort described in this report has been directed at developing 
fundamental analytical models of the heave and the heave- pitch motion of 
Air Cushion Landing Gear. These models have been implemented in a 
computer program delivered to NASA. The program is capable of 
simulating the dynamic heave and heave -pitch behavior (aircraft g load- 
ing and motion, trunk deflection, pressures, etc. ) of an AC LS- equipped 
aircraft caused by landing impact and taxi over an irregular runway, 
using input data such as aircraft weight, ACLS geometry, fan character- 
istics, runway surface profile, etc. Three types of illustrative simula- 
tions have been carried out to demonstrate the capabilities of the pro- 
gram. The illustrative results show how drop tests, landing impact 
and rough runway operation can be simulated. 

In the next phase of this program, a coupled heave -pitch- roll 
analysis will be developed. Also, experimental verification and refine- 
ment of the analysis using test data obtained at NASA- Langley with a 
model cushion will be performed. After the program capabilities have 
been verified, more extensive simulations are planned, to investigate 
a variety of potentially attractive cushion configurations and to develop 
guidelines for improved ACLS designs. 
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APPENDIX A - PROGRAM ORGANIZATION AND USE 


The overall structure of the computer program developed for 
simulating the heave-pitch dynamics of the air cushion landing systems 
is described in this Appendix, along with instructions on its usage. 
Appendices B, D, E and F described various aspects of the program 
in greater details. 

A. 1 Program Organization 

The ACLS heave -pitch model is simulated as follows; 

(a) The input data is read in. 

(b) Initial geometry calculations are carried out. 

(c) The static characteristics are computed and printed. 

(d) The initial conditions for the dynamic simulation are 
determined. 

(e) The state equations are integrated numerically to 
determine the time history of ACLS pressure and 
motion following landing impact. 

The computer program that simulates ACLS heave-pitch dynamics 
is listed in Appendix E. The flow diagram is shown in Figure A. 1, 
and the computing sequence is described below. 

(a) Data Input and Conversion 

Initially, subroutine PROGIO is called which reads the 
input data cards through five other subroutines. Each card contains 
alphanumeric data which includes the name of each parameter, its value, 
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and dimension. The input parameters are printed directly after they 
are read. Some of the less frequently altered parameters are specified 
directly in the main program. The input data is then converted, in this 
case, to ft-lb-sec units prior to the computations. 

(b) Initial Geometry Calculations 

Subroutine TRUNK calculates the trunk shape parameters 
(radii of curvature, subtended angles, etc. ) from the input parameters 
i, hy, a, b, d and L g (see Figure 1). Subroutine SEGMNT divides the 
trunk into a (user specified) number of segments and calculates the seg- 
ment center distance from the center of the cushion. Subroutine SHAPE 1 
then calculates the trunk cross section area, trunk volume, cushion area, 
trunk-to- cushion orifice area, trunk-to-atmosphere orifice area, and 
distance of the center of cushion pressure for each segment, when the 
trunk is out of ground contact. 

These three subroutines, TRUNK, SEGMENT and SHAPE 1, 
are called only once in the program. The values of areas and volumes 
assessed by SHAPE 1 are independent of ACLS motion. Another sub- 
routine, SHAPEZ, is called in the program whenever updated values of 
areas and volumes are required. 

(c) Static Characteristics 


Subroutine FLOW1 is called next. This subroutine cal- 
culates the static characteristics, i. e. , the height of the aircraft center 
of gravity from the ground, pitch angle, gap area, plenum, trunk and 
cushion pressures, and total flow for various combinations of aircraft 
load and position of the center of gravity. This is accomplished by 
calling six subroutines; COORDN, PROFILE, CLRNCE, SHAPE2, FORCE 
and CMPCRV. The first four subroutines calculate the required areas 
and volumes of the ACLS for a particular combination of CG height and 
pitch angle. Subroutine FORCE calculates the torque and load developed 
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Figure A. 1 Program Flow Diagram - Initialization 
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by the particular configuration and CMPCRV supplies the fan character- 
istics to F LOW 1 . 

The static characteristics are used to determine the final 
equilibrium values of the aircraft CG height above the ground, pitch 
angle, areas, volumes, pressures, flows and stiffness i$or the input 
value of aircraft weight and CG position. The static characteristics 
(10 values) and the final equilibrium values are then printed. 

(d) Initial Conditions 


Initial values of the state variables are needed to start 
the Runge-Kutta integration in the dynamic simulation. The initial 
values of CG coordinates, pitch angle, sink rate, horizontal velocity 
and pitch (rotational) velocity are supplied by the user. Subroutines 
COORDN, PROFILE, CLRNCE and SHAPE2 are called to determine the 
initial cushion and trunk geometry. Then the initial values of the 
cushion pressure, trunk pressure and plenum pressure are calculated 
from the static characteristics, computed as described in (c) above. 

(e) Integration of State Equations 

The dynamic simulation is coordinated by subroutine 
DYSYS. DYSYS calls subroutine RKDIF at each time step - RKDIF 
starts with the values of the seven state variables (plenum, trunk and 
cushion pressure, CG height, sink rate, pitch angle and pitch velocity) 
at a given time (t), and calculates new values of these variables at 
time (t+dt) by numerical integration using the Runge-Kutta method. The 
new values are obtained from derivations of the state variables at time 
t, and at intermediate times between t and t+dt. The calculation of 
derivatives is carried out in subroutine STEQU. STEQU contains the 
basic differential equations of pressure and motion of the ACLS. In 
order to calculate the derivatives, STEQU requires values of flows, 
forces and torques fbr the given values of the state variables. This is 


51 



accomplished by calling subroutines COORDN, PROFILE, CLRNCE, 
SHAPE2, FLOW2 and FORCE, in that order. The first four subroutines 
calculate the various trunk and cushion areas and volumes correspond- 
ing to the ACLS orientation at a particular instant of time. From this 
data, subroutine FLOW2 determines the various flows through the ACLS. 
Subroutine CMPCRV supplies FLOW2 with data on the pressure-flow fan 
characteristics. Finally, subroutine FORCE calculates the forces and 
torques acting on the aircraft eg using the values of the appropriate 
pressures and areas supplied by STEQU and SHAPE2 respectively. 

The above procedure is repeated for each time increment, 
and output values of pressures, flows, motion, etc. , are printed at 
appropriate (user specified) intervals. The program terminates when 
the simulation time equals the user supplied time limit. 
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A. 2 Program Use 


A. 2. 1 Input Data Format 

Input data is supplied to the program in three ways. 

(a) By data cards that are read in after program 
compilation (i. e. , through the READ statement). 

(b) By data specifications that are included within 
the program (i. e. , through DATA and other 
specification statements ). 

(c) Through Subroutine PROFILE which is used to 
specify ground profile information. 

Method (a) is used primarily to specify the parameters 
that are design variables and/or are likely to be frequently changed 
(aircraft weight, initial sink rate, initial pitch angle, etc. ). Method (b) 
is used primarily for parameters which are likely to be changed less 
frequently (discharge coefficients, polytropic exponent, etc. ). The for- 
mat for specifying these two types of input data are given below. The 
input data are in English units. 

The profile description essentially consists of supplying 

values of ground elevation Y (i) corresponding to segment i (Eq. (19), 

8 

Appendix C). This can be accomplished in one of several ways, e. g. , 
using function subprograms, data cards, etc. , depending on user prefer- 
ence and form in which the data is available. 

A. 2. 1. 1 Data Cards 

The following data cards are required to execute 
the program. Most cards have an alphanumeric input format. The 
numerical values are placed within the T F ', ■I 1 or 'E 1 fields. The ! A' 

fields are used to provide parameter names and units, and other legends 
for user convenience. 
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Card No. 

Contents 


Format 

1-7 

Header cards 

80A1 




(The user can print a heading using 
these cards) 





Aircraft Parameters 




8 

Aircraft weight (lbs) 

30A1, 

F10. 4, 

10A1 

9 

■ Aircraft pitch inertia about CG 
(slug ft 2 ) 

30A1, 

F10. 4, 

10A1 

10 

Horizontal distance of CG from 
geometric center of cushion, CC (ft) 

30A1, 

F10. 4, 

10A1 

11 

Vertical distance of CG from 
geometric center of cushion, GG (ft) 

30A1, 

F 10. 4, 

10A1 

12 

Heave drag coefficient, C^ 

30A1, 

F10. 4, 

10A1 

13 

Heave drag area, A ^ (ft^) 

30A1, 

F 10. 4, 

10A1 

' 

Trunk Parameters 




14 

Header card 

80A1 



15 

Straight section length of cushion, 

L s (ft) 

30A1, 

FI 0. 4, 

10A1 

16 

Distance between inner trunk 
attachment points, d (ft) 

30A1, 

F 1 0. 4, 

10A1 

17 

Horizontal distance between inner 
and outer trunk attachment points. 

30A1, 

F10. 4, 

10A1 


. a (ft) 




18 

Vertical distance between inner and 
outer trunk attachment points, b (ft) 

30A 1, 

F10. 4, 

10A1 

19 

Peripheral length of trunk cross 
section, i (ft) 

30A1, 

F 10. 4, 

10A1 

20 

Inflated trunk height (no load), 
h (ft) 

30A1, 

F10. 4, 

10A1 


y 




2 1 


"3 A A 1 

T C 


22 

' - — - - * r 

Number of trunk orifices per row, N^ 

30A1, 

15 
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Format 


Card No. 

23 

24 

25 

26 

27-28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


Contents 


2 

Area of each trunk orifice, (in ) 
Spacing between trunk orifice rows, 

s h (ft) 

Header card 


30A1, 

30 . 1 . 1 , 


F10. 4, 
F10. 4, 


10A1 

10A1 


80A1 


Peripheral distance between inner trunk 
attachment point and first row of holes. 


L (ft) 
P 

30A1. F10. 4, 

10A1 

Air Supply Parameters 



Header cards 

80A1 


Plenum-to- cushion orifice area, 
A plch (ft2 > 

30A1, F10. 4, 

10A1 

Plenum-to-trunk orifice area, 
A P itk < ft2 > 

30A1, F10. 4, 

10A1 

Plenum-to- atmosphere orifice area, 
A plat < ft2 > 

30A1, F10. 4, 

10A1 

2 

Fan inlet orifice area, A , (ft ) 
(See Note 1) atln 

30A1, F10. 4, 

10A1 

3 

Plenum volume, (ft ) 

30A1, F10. 4, 

10A1 

3 

Dead cushion volume, V ^ (ft ) 

30A1, F 10. 4, 

10A1 

Landing Approach Parameters 



Header card 

80A1 


Initial X coordinate of CG, XCGI (ft) 

30A 1, F10. 4, 

10A1 

Initial Y coordinate of CG, YCGI (ft) 

30A1, F10. 4, 

10A1 

Initial pitch angle, PHII (degrees) 

30A1, F10. 4, 

10A1 

Initial horizontal velocity, VELXI 
(ft/sec) 

30A1, F10. 4, 

10A1 

Initial sink rate, SINKRT (ft/sec) 

30A1, F10. 4, 

10A1 

Initial pitch velocity, DPHI 
(degrees /sec) 

30A1, F10. 4, 

10A1 
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Card No. Contents Format 

Environmental Conditions 


42 

Header card 


80A1 



43 

Absolute atmospheric pressure, 

P at (psia) 

30A1, 

F10. 4, 

10A1 

44 

Ambient temperature (°F) 


30A1, 

F10. 4, 

10A1 


Fan Characteristics 





45 

Fan curve polynomial coefficient* & Q 

40A1, 

E10. 3 


46 

II 

a l 

40A1, 

E10. 3 


47 

II 

a z 

40A1, 

E10. 3 


48 

II 

a 3 

40A1, 

E10. 3 


49 

It 

a 4 

40A1, 

E10. 3 


50 

M 

Po 

40A1, 

E10. 3 


51 

M 

01 

40A1, 

E10. 3 


52 

II 

02 

40A1, 

E10. 3 


53 

II 

03 

40A1, 

E10. 3 


54 

II 

04 

40A1, 

E10. 3 


55 

Maximum unstalled fan pressure rise, QP1 

40A1, 

E10. 3 


56 

Minimum stalled fan pressure 

rise, QP2 

40A1, 

E10. 3 


57 

Minimum unstalled flow, QP3 


40A1, 

E10. 3 


58 

Maximum unstalled flow, QP4 


40A1, 

E10. 3 


59 

4 

Maximum stalled flow, QP5 


40A1, 

E10. 3 



Integration and Print Control Parameters 




60 

Integration time step, DTIME 

(sec) 

40A1, 

E10. 3 


61 

Simulation time limit, FTIME 

(sec) 

40A1, 

E10. 3 


62 

Number of time steps between 
printing, MM 


40A1, 

15 


See 

Figure 7 for explanation of symbols 

on cards 45 

-59. 
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Card No. 


Contents 


Format 


Trunk Segment Parameters 

63 Number of straight sections in one-fourth 

of the periphery, M 40A1, 15 

64 Number of curved sections in one-fourth 

of the periphery, N 40A1, 15 

Note 1. Care must be taken to choose the proper area for the fan 

upstream orifice. When the orifice is far enough upstream so that it is 
not affected by the fan inlet flow field, the effective orifice area can be 
found from geometry. When the upstream orifice is close to the fan inlet 
(e. g. , partial inlet blockage), the flow patterns in the fan inlet will affect 
the orifice characteristics, and the effective orifice area must be found 
from measurements of the flow and pressure drop. For values of A^ tfa 
larger than 1 ft 2 , the simulation is carried out with an unrestricted fan 
inlet (no upstream orifice). The value 1 ft 2 is chosen arbitrarily and it can be 
altered, if necessary, by changing Statement No. 76 in Subroutine FL.OW2 on 
page 177. 

A. 2. 1. 2 Internal Data 

The internal data are specified at the beginning 
of the main program (see line marked 'DATA ACQUISITION', Page 149, 


Appendix E). They are 

as follows. 


CPA 

Discharge coefficient, plenum-to- 
atmosphere orifice 

CAF 

Discharge 

coefficient, fan inlet orifice 

CPC 

Discharge 

orifice 

coefficient, plenum-to- cushion 

CPT 

Discharge 

orifice 

coefficient, plenum- to -trunk 

CTC 

Discharge 

orifice 

coefficient, trunk- to- cushion 

CGAP 

Discharge 

coefficient of clearance gap 
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CTA 


Discharge coefficient, trunk- to- atmosphere 
orifice 


CKK 

- 

Polytropic expansion exponent 

GEC 

- 

Ground effect coefficient (See Note 1) 

ZETA 

- 

Trunk damping ratio (See Note 2) 

PERTK 

- 

Trunk orifice blockage parameter (See 
Note 3) 

PERCH 

- 

Cushion exit seal parameter (See Note 3) 

U 


Ground-trunk friction coefficient (See 
Note 4) 

DECCL 

- 

/ 2 

Aircraft horizontal deceleration rate (ft/sec ) 
(See Note 5) 


Note 1. When the ACLS is high above the ground, the cushion 

pressure (gage) will be zero. Although simulation of the full dynamic 
model will predict this condition, significant computing economy can be 
achieved by simplifying the full model at large heights by assuming zero 
cushion pressure rather than obtaining this same result through solution 
of the differential equation of cushion pressure. The ground effect coeffi- 
cient is the factor which determines the ground effect area (cushion-to- 
atmosphere gap) above which the cushion pressure is set equal to zero 
rather than computed from the full ACLS simulation. The ground effect 

gap area, A , is determined from the ground effect coefficient as 
s gapg 

follows 


gapg 


(GEC) 


where is the equilibrium gap area found from the static load- 

deflection characteristic. Initial simulations with GEC = 10 have given 
satisfactory results. Larger values will require smaller integration time 
steps, and involve more computation. Smaller values may allow larger 


58 


time steps, but can lead to starting transients when the AC LS comes 
into ground effect. 


Note 2. The trunk damping ratio zeta, f , is a nondime ns ion al measure 

of trunk damping. The damping coefficient for each segment, B , 
[Assumption (f) in Section Z.z\ is obtained from the damping ratio as 
follows 


B = Zt \/ Heave Stiffness x Mass /4(M+N) 
z 

This assumes that the damping force is equally divided amongst all trunk 
segments. Test data and dimensional analysis will provide the basis for 
estimating the trunk damping ratio. Initial data indicates that the damp- 
ing ratio will be in the range of 0. 05 - 0.1. 

Note 3. During ground contact, the trunk orifices nominally covered 

by the ground will not be completely blocked, and the cushion will not be 
perfectly sealed, because of ground irregularities, trunk ribs and imper- 
fections and brake pads. Therefore, some small flow will leak out the 
cushion to the atmosphere and through the covered trunk orifices even 
during ground contact. PERTK and PERCH are measures of the trunk 
and cushion leakage areas during ground contact. PERTK is the fraction 
of the trunk orifice area that is blocked during ground contact. For 
example, PERTK = . 85 signifies that 85% of the area of the trunk orifices 
in ground contact is blocked, and the leakage flow occurs through the 
unblocked 15% of the orifice area. PERCH is the ratio of the cushion 
leakage area during ground contact to the equilibrium clearance gap area. 
For example, PERCH = . 15 signifies that during ground contact, 15% of 
the equilibrium clearance gap area remains unblocked. 

Typically, when the brake pads are not deployed, PERTK 
will be aDout .85 - .9 and PERCH will be about .1 - .15 (i. e. , 85-90% 
blockage in both cases). When the brake pads are deployed, the blockage 
will be reduced, and appropriate values for PERTK and PERCH can be 
estimated from pad geometry. 
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Note 4. The trunk- ground friction coefficient is required to 

calculate the pitching torque due to friction force. It depends on many 
factors, including the brake pads and/or trunk material characteristics, 
ground surface characteristics, etc. For simulations carried out thus 
far, the friction coefficient has been assumed to lie in the range of 

0. 4 to 0. 5. 


Note 5. In this simulation, the aircraft is assumed to decelerate, 

in a horizontal direction, at a constant (user selected) rate. 

A. 2. 2 Program Output 


dix F. ) 


The output printout provides the following data. (See Appen- 

1. The Input Parameters 


Aircraft weight and pitch inertia, location of eg, 
initial approach parameters (sink rate, etc. ), ACLS configuration, etc. 


2. Final Equilibrium Conditions 


(These are the conditions that exist after the 
landing dynamics have damped out. ) 

Height of CG 
Pitch Angle 
Cushion Perimeter 
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Cushion Volume 
Trunk Volume 

Gap Area, Cushion- atmosphere 
Cushion Area 
Trunk Contact Area 
Orifice Area, Trunk-atmosphere 
Orifice Area, Trunk-cushion 
Cushion Pressure (gage) 

Trunk Pressure (gage) 

Plenum Pressure (gage) 

Total Volume Air Flow 
Total Volume Cushion Flow 
Volume Flow, Plenum to Cushion 
Volume Flow, Plenum to Trunk 
Volume Flow, Trunk to Cushion 
Volume Flow, Trunk to Atmosphere 
Volume Flow, Plenum to Atmosphere 
Stall Margin (See Note 1) 

Heave Stiffness 
Pitch Stiffness 
Theoretical Fan Power 

3. Static Characteristics 


The height of the center of gravity above 
the ground, pitch angle, air gap area, cushion and trunk pressures and 
total flow for various combinations of force and torque loads (i. e. , for 
various weights and offset distances from the center of the cushion). 

4. Dynamic Simulation 

Evaluation of the following variables at successive 
intervals of time during aircraft approach, touchdown and taxi. 
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Acceleration (vertical) 

Velocity (sink rate) 

CG Height (Y coordinate of eg) (from reference 
X axis) 

X coordinate of CG 
Trunk Pressure (gage) 

Cushion Pressure (gage) 

Fan Volume Flow 
Pitch Angle 

Trunk- to -Cushion Volume Flow 

Trunk-to-Atmosphere Volume Flow 

Cushion- to- Atmosphere Volume Flow 

Gap Area (clearance area between trunk and ground) 

Note 1. The fan stall margin is the maximum percentage rise in fan 

pressure that can occur without fan stall. This parameter is a measure 
of the ability of the system to absorb dynamic impact without fan stall. 

When the stall margin is below 5% (i. e. , impact stall likely), the state- 
ment 'FAN CRITICALLY STABLE' appears in the program output. 

A. 2. 3 Premature Program Termination 

Premature program termination can occur under the follow- 
ing conditions: 

(a) When the input parameters do not allow feasible 
solutions to be obtained. For instance, when 
the aircraft weight cannot be supported due to 
insufficient output from the fan. 

(b) When the integration time increment DTIME is 
chosen too large, and causes numerical 
instabilities during the solution of the differential 
equations. 
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A. 2. 3. 1 Diagnostic Messages 


When input parameter values do not allow a 
feasible solution, one of the following diagnostic statements is printed 
prior to program termination. 

1. 'INFEASIBLE TRUNK GEOMETRY' 

This message indicates that the geometrical 
trunk parameters a, b, f and h^ are such 
that a feasible solution for the trunk shape 
cannot be obtained - i. e. , a real curve 
joining the trunk attachment points, and 
having length i and height h^ cannot be 
found. 

2. 'INFEASIBLE CONFIGURATION' 

This message indicates that the ACLS con- 
figuration specified by the user is infeasible. 
The means to correct this situation include; 

(a) Reduced load 

(b) Increased fan output 

(c) Reduced plenum bleed area 

3. 'CHANGE VALUE OF XTOL' 

In the rare case that this message appears, 
the situation is remedied by increasing the 
iteration tolerance value XTOL. 
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A. 2. 3. 2 Numerical Instability 


When the user supplied integration time step 
DTIME is too large, the Runge Kutta integration scheme will not con- 
verge, and the program will terminate due to field overflow. In such 
cases, a smaller time step will eliminate the problem. Initial results 

indicate that for typical ACLS configurations and operating conditions, 

-4 

time steps less than about 5 x 10" sec give satisfactory results. 

A. 2. 4 Heave Simulation 


The heave-pitch simulation program can be adapted for 

simulating just the heave motion by assuming a very high value for 

10 2 

aircraft pitch inertia (say, 10 slug ft ). The pitch inertia data are 
supplied by card No. 9, as described in Section A. 2. 1. 1. Initial 
pitch angle and pitch velocity (Card No. 38 and 41, respectively) should 
be zero, while simulating just the heave motion. 
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APPENDIX B - PRINCIPAL PROGRAM NOMENCLATURE* 


The symbols used in the analysis and the corresponding computer 
program variable names are defined below. All program variables are 
computed in the appropriate ft-lb-sec units except where indicated to the 


contrary. 


• 

Symbol 

Program 
Variable Name 

Explanation 

a 

A 

Horizontal distance between inner 
and outer trunk attachment points 

A r A 9 

A1-A9 

Area of trunk sectors used in 
volume computations 

A atfn 

AATFN 

Fan inlet orifice area 

- 

ACCEL 

Aircraft acceleration (positive upwards) 

A ch 

ACH 

Cushion area 

A chi<‘> 

ACHI(I) 

I value of cushion area (i. e. , with 
no ground contact) of ith segment 

A chr<‘> 

ACHR(I) 

R value of cushion area (i. e. , 
decrement due to ground contact) 
of ith segment 


AD IF 

Gap area above ground effect gap 
area 

A gap 

AGAP 

Clearance gap area 

A gape 

AGAPE 

Equilibrium gap area 

A gapi (l) 

AGAPI(I) 

I value of clearance gap area (i. e. , 
with no ground contact) of ith segment 

- 

AGAPP(I) 

Iterated gap area 

A g>pr (i » 

AGAPR(I) 

R value of clearance gap area (i. e. , 
decrement due to ground contact) 
of ith segment 

- 

AG APS (I) 

Static values of AGAP 


AH 

Area of trunk orifice 


^'Nomenclature is listed alphabetically by symbols. Therefore some 
program variable names do not appear alphabetically. 
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Symbol 

Program 
Variable Name 

Explanation 

A leak 

ALEAK 

Minimum clearance gap area 
(caused by brake pads, imperfections, 
etc. ) 

A ph 

PHA 

Heave drag area of cushion 

A plat 

APLAT 

Plenum -to -atmosphere orifice area 

A plch 

APLCH 

Plenum- to -cushion orifice area 

A pltk 

APLTK 

Plenum-to-trunk orifice area 

A tk 

ATK 

Trunk cross sectional area 

A tki (i) 

ATKI(I) 

I value of trunk cross sectional area 
of ith segment 

A tkr<*> 

ATKR(I) 

R value of trunk cross sectional area 
of ith segment 

A tkat 

ATKAT 

Trunk- to -atmosphere orifice area 

^kati^ 

ATKATI(I) 

I value of trunk-to-atmosphere 
orifice area of ith segment 

A tkatr<^ 

ATKATR(I) 

R value of trunk-to-atmosphere 
orifice area of ith segment 

A tkch 

ATKCH 

Trunk-to-cushion orifice area 

Atkchi^^ 

ATKCHI(I) 

I value of trunk-to-cushion orifice 
area of ith segment 

Atkchr^) 

ATKCHR(I) 

R value of trunk-to-cushion orifice 
area of ith segment 

A., 

tkcn 

ATKCN 

Trunk- ground contact area 

A tkcni^ 

ATKCNI(I) 

I value of trunk-ground contact 
area of ith segment 

A tkcnr^ l) 

ATKCNR(I) 

R value of trunk- ground contact 
area of ith segment 

- 

ATOL 

Tolerance for gap area iteration 

b 

B 

Vertical distance between trunk 
attachment points 
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Symbol 

Program 
Variable Name 

Explanation 

B 

z 

DAMPC 

Trunk damping constant 

B 

z 

- 

Damping constant for each trunk 
segment 

C af 

CAF 

Discharge coefficient of atmosphere - 
to-fan orifice 

CC 

CC 

Horizontal distance of CG from 
center of cushion 

- 

CCI 

User supplied value of CC 

' - 

CCS(I) 

Values of CC used in the static 
computations 

C D 

HDC 

Heave drag coefficient 

C enf 

CENF 

Distance of center of aerodynamic 
heave drag area from CG 

C 

gap 

CGAP 

Discharge coefficient of clearance 
gap 

C 

pa 

CPA 

Discharge coefficient of plenum-to- 
atmosphere orifice 

C pc 

CPC 

Discharge coefficient of plenum-to- 
cushion orifice 

C pt 

CPT 

Discharge coefficient of plenum-to- 
trunk orifice 

C ta 

CTA 

Discharge coefficient of trunk-to- 
atmosphere orifice 

C tc 

CTC 

Discharge coefficient of trunk-to- 
cushion orifice 

d 

D 

Distance between trunk attachment 
points 

- 

DECCL 

Braking deceleration 

del 

X 

DELX 

Width of straight trunk segment 
(along aircraft axis) 

. - 

DER(I) 

Derivatives of state variables 
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Symbol 


Program 
Variable Name 


Explanation 


• 

DPHI 

Pitch (angular) velocity 

- 

DPHII 

Initial pitch velocity 

- 

DTIME 

Time step dt 

- 

DVCH 

Time derivative of cushion volume 

- 

DVTK 

Time derivative of trunk volume 

- 

DY(I) 

Derivatives of state variables 

- 

FORCD 

Aerodynamic drag force 

Forcn 

FORCN 

Total vertical pressure force 
transmitted to aircraft 

- 

FORCNS(I) 

Static values of FORCN 

Foret 

FORCT 

Trunk damping force 

- 

FTIME 

Simulation time limit 

- 

GEC 

Ground effect coefficient 

GG 

GG 

Vertical distance of CG from center 
of cushion 

- 

HP 

Theoretical fan power 

h 

y 

HY 

Equilibrium height of trunk cross 
section 

- 

ICASE 

Element number of iteration closest 
to the given gap area. ICASE = 0 
if configuration infeasible. 

- 

ICLN 

Used to suppress error message of 
CLRNCE during iteration of FLOW1 

- 

ICON 

Controls storage of static character- 
istic values in FLOW1. 

- 

ICREST 

ICREST = 2 if configuration load and 
position within tolerance to aircraft 
weight and CG position. 

ICREST = 1 otherwise 
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Program 
Variable Name 


Explanation 


ID 

ID ID 

IDIF 
IF AN 

HT 

INERT 

IODIN 

IPCT 

IPHI 

IPN 

IPP 

IPREST 

IQ 

IRST 

IS 


ID = 0 after first entry in ground 
effect region. ID = 1, cushion out 
of ground effect (initially) 

IDID = 0 if any configuration of 
iteration 1 in FLOW1 infeasible. 

IDID = 1 otherwise 

Same as ICASE 

IF AN = 0 for unstalled fan operation 
IF AN = 1 for stalled fan operation 

Element number in static character- 
istic 

Pitch moment of inertia of aircraft 
about CG 

Indicator of number of passes 
through iteration 1 of FLOW1 

IPCT = 0 if PCH > PTK 
IPCT = 1 if PCH< PTK 

Y element number of two dimensional 
array of iteration 1 grid in FLOW1 

Indicator of number of passes 
through iteration 2 of FLOW1 

IPP = 0 for combined trunk, plenum 
dynamics, IPP = 1 separate trunk, 
plenum dynamics 

IPREST = 2 if actual gap area and 
iterated gap area within tolerance 
IPREST = 1 if parameters not 
within tolerance 

Used to effect gradual change in 
cushion flow model after ground 
effect zone entry 

Counter used to prevent endless 
iteration in FLOW1 

IS = IXCG of the best iteration 
grid point 


Symbol 

Program 
Variable Name 

Explanation 

- 

ISAVE 

Same as ICASE 

- 

ISHAPE 

ISHAPE = 0 if trunk inflation impossible 
ISHAPE = 1 if trunk inflation possible 

- 

IYCG 

X element number of two dimensional 
array of iteration 1 grid in FLOW1 

- 

JS 

JS - IPHI of the best iteration grid ^ 

point 

a 

L 

Peripheral trunk length 

h 

LI 

Trunk length, outer attachment to 
horizontal tangent 

h 

L2 

Trunk length, inner attachment to 
horizontal tangent 

i 

p 

LP 

Peripheral distance from inner 
trunk attachment to first row of 
orifices 

L 

s 

LS 

Straight section length of cushion 

M 

M 

Number of straight trunk segments 
in one quarter of trunk periphery 

M 

a 

MASS 

Mass supported by ACLS (slugs) 

- 

MM 

Number of steps before printing 

N 

N 

Number of curved trunk segments 
in one quarter of trunk periphery 

N h 

NH 

Number of trunk orifices per row 

- 

NQ 

Determines number of steps in 
which gap model change takes place 

N 

r 

NR 

Number of rows of trunk orifices 

P at 

PAT 

Atmospheric pressure, absolute 

P atfn 

PATFN 

Fan inlet pressure, gage 

P ch 

PCH 

Cushion pressure, gage 
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Program 



Variable Name 

Explanation 

- 

PCHS(I) 

Static values of PCH, gage 

- 

PEN(J) 

Iteration value of PATFN 

Per ch 

PERCH 

Determines air gap seal imperfection 

Per tk 

PERTK 

Factor that determines flow area of 
trunk orifices in ground contact 

P fan 

PFAN 

Fan pressure rise 

“ 

PFANS(J) 

Iterative values of PFAN 

- 

PHIDELT 

Increment in pitch angle 

" 

PHH 

Initial value of pitch angle 


PHIS (I) 

Static values of pitch angle 

“ 

PHISTRT 

Initial value of pitch angle in 
iteration 1 of FLOW1 

- 

PHISTOP 

Final value of pitch angle in 
iteration 1 of FLOW1 


PINC 

Increment in PFAN after first 
feasible configuration 

“ 

PINC1 

Increment in PFAN before first 
feasible configuration 

P plm 

PPLM 

Plenum pressure, gage 

“ 

PPLMS(I) 

Static values of PPLM, gage 

- 

PSTART 

Starting value of PFAN for iteration 

- 

PS TOP 

Last value of PFAN in iteration 

P tk 

PTK 

Trunk pressure, gage 

- 

PTKS(I) 

Static values of PTK, gage 

- 

PTOL 

Tolerance for pressure and force 
iteration 

^chat 

QCHAT 

Cushion-to-atmosphere flow (cfs) 
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Program 


Symbol 

Variable Name 

Q fan 

QFAN 

- 

QFANS(I) 

Q fnpl 

QFNPL 

- 

QP1-QP5 

Q plat 

QPLAT 

^plch 

QPLCH 

Q pltk 

QPLTK 

^tkat 

QTKAT 

^tkch 

QTKCH 

- 

RTOL 

R 1 

R1 

R 2 

R2 

S 

S 

- 

SC 

S h 

SH 

- 

SINKRT 

- 

SINPH2 

- 

SINPHR 

T 

cp 

- 

T 

df 

- 

- 

TEMP AT 

t 

TIME 

T. 

tp 

- 



Fan flow (cfs) 

Static values of QFAN 

Fan-to-plenum flow (cfs) 

Stall, recovery and other points on 
fan characteristic 

Bleed flow (cfs) 

Plenum- to -cushion flow (cfs) 

Plenum -to -trunk flow (cfs) 

Trunk- to -atmosphere flow (cfs) 

Trunk-to-cushion flow (cfs) 

Tolerance for R 2 iteration 

Outer radius of curvature of trunk 

Inner radius of curvature of trunk 

Peripheral length of cushion 

Stall margin 

Trunk orifice row spacing 
Aircraft velocity (positive upward) 
Sin 

R 2 Sin 0 2 

Cushion pressure torque component 
Torque due to drag force 
Ambient temperature 
Time 

Trunk pressure torque component 
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Program 


Symbol 

Variable Name 

Explanation 

Torf 

TORF 

Trunk- ground contact friction torque 

Torn 

TORN 

Pressure torque 

- 

TORQ 

= TORN + TORF 

Torqt 

TORQT 

Trunk damping torque 

V ch 

VCH 

Total cushion volume 

V chd 

VCHD 

Cushion dead volume 

V chi< i » 

VCHI(I) 

I value of (cushion volume of ith 
segment 

V chr<‘» 

VCHR(I) 

R value of cushion volume of ith 
segment 


VCHS 

VCH at time t-dt 


velx 

Forward velocity of the aircraft 

- 

VELXI 

Initial forward velocity 

V plm 

VPLiM 

Plenum volume 

V tk 

VTK 

Trunk volume 


VTKI(I) 

I value of trunk volume of ith 
segment 


VTKR(I) 

R value of trunk volume of ith 
segment 

- 

VTKS 

VTK at time t-dt 

X 

X 

Variable used in trunk area calculations 

x r x 9 

X1-X9 

X-coordinates of centroids of trunk 
cross sectional area components 

X 12 

X12 

X coordinate of trunk area components 
(A1-A2) 

X cc 

XCC 

X- coordinate of center of cushion (CC) 

X cg 

XCG 

X-coordinate of CG 
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Program 


Symbol 

Variable Name 

Explanation 

- 

XCGI 

Initial XCG 

X ch< i > 

XCH(I) 

Distance of center of cushion 
pressure of ith segment from CC 

X cr 

XCR 

X- coordinate of trunk area components 
(A6-A7) 


XCX(I) 

Distance of center of ith segment 
from CC 

X e 

XE 

X- coordinate of centroid of total 
trunk cross section area 

X or 

XER 

X-coordinate of centroid of area 
A tkr (l) 


XH(I) 

X-coordinate of center of ith 
segment 

X tk (i) 

XTK(I) 

Distance of center of trunk pressure 
of ith segment from CC 

- 

XTOL 

Tolerance for iteration 1 of FLOW1 

- 

Y(I) 

State variables (see below) 

- 

Y(l) 

Plenum pressure, gage 

- 

Y(2) 

Cushion pressure, gage 


Y(3) 

Trunk pressure, gage 

- 

Y(4) 

Vertical aircraft velocity 

- 

Y(5) 

Vertical aircraft displacement 

- 

Y(6) 

Pitch (angular) velocity 

• 

Y(7) 

Pitch angle 

Y cc 

YCC 

Y-coordinate of center of cushion 
(CG) 

Y cg 

YCG 

Y-coordinate of CG 

- 

YCGDELT 

Increment in YCG 
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Program 

Symbol Variable Name Explanation 


- 

YCGI 

Initial value of YCG 

- 

YCGS(I) 

Static values of YCG 

- 

YCGSTOP 

Final value of YCG in iteration 1 
of FLOW 1 

- 

YCGSTRT 

Initial value of YCG in iteration 1 
of FLOW 1 

- 

YD IF 

= ADIF/S 

Y,(i) 

YG(I) 

Ground elevation corresponding to 
ith segment 

V‘> 

YGH(I) 

Hard surface clearance for ith 
segment 

Y h<‘) 

YH(I) 

Y-coordinate of center of ith 
segment 


ZCC(I, J) 

Value of CC for point (I, J) in ' 

iteration 1 of FLOW1 

- 

ZWT(I, J) 

Value of FORCN for point (I, J) 
in iteration 1 of FLOW1 

“O' “4 

AL0-AL4 

Polynomial coefficients of fan curve 

0 

BETA 

Angle subtended by curved segment 
of the trunk 

"CD 

0 

1 

Xo 

4* 

B0-B4 

Polynomial coefficients of fan curve 

a (i) 

DELTA(I) 

Angular position of ith curved 
segment 

K 

CKK 

Polytropic expansion constant 

0 

PHI 

Pitch angle. Positive anticlockwise 

0 1 

PHI1 

Angle subtended by outer trunk 
segment (atmosphere side) 

0 2 

PHI2 

Angle subtended by inner trunk 
segment (cushion side) 

0 3 

PH 13 

Angle subtended by cushion side of 
trunk deformation (Figure 6) 
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Symbol 


P 

P 

t 


Program 
Variable Name 

Explanation 

PHI4 

Angle subtended by atmosphere side 
of trunk deformation (Figure 6 ) 

RHO 

3 

Air density (slugs /ft ) 

U 

Coefficient of friction between the 
trunk and the ground 

ZETA 

Trunk damping ratio 
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APPENDIX C - DETAILED HEAVE-PITCH MODEL ANALYSIS 


C. 1 Introduction 


The heave-pitch model incorporated in the program is divided into 
two parts - the static relationships and the dynamic model. 

. The static relationships evaluate the geometric parameters, pres- 
sures and flows for the ACLS in equilibrium. They provide static design 
data, and also determine the initial conditions for the dynamic simulation. 

The dynamic model predicts the time histories of ACLS pressures 
and flows and aircraft motion during approach, touchdown and taxi. 

C . 2 The Static Model 


The static relations . can be divided into two categories; 

(a) Geometric relations, summarized in C. 2. 1 and C.2. 2. 

(b) Pressure-flow-force-torque relations, summarized in C.2. 3 

C. 2. 1 Trunk Crossectional Shape 

From trunk geometry (Figure 6), 


R 1 0 1 =i l 


R 2*2 =i 2 


*1 + h = * 


R_-h 

Cos *2 = 


( 1 ) 

( 2 ) 

(3) 


(4) 


77 


a 


(5) 


RjCos^-gO 0 ) + R 2 Sin0 2 = 

RjSin^-W 0 ) - (h -Rj) = 1? (6) 

The six unknown trunk configuration parameters ly 
Rj, R^, and 0^ can -be obtained by solving Equations (1) through (6) 
simultaneously, in terms of the known trunk parameters a, b, h^ and L 

C. 2. 2 Segmented Trunk Analysis 

The orifice areas and cushion and trunk volumes, for a 
particular trunk orientation, are calculated as follows; 


The trunk is divided into a number of segments. Then 
for a given trunk orientation, areas and volumes are calculated independ- 
ently for each segment. The areas and volumes are divided into two com- 
ponents - the i component and the r component. The i values are calcu- 
lated assuming that the trunk segment under consideration is out of ground 
contact. The r values represent the changes of the segment areas and 
volumes due to trunk-ground contact. The segment areas and volumes are 
found by subtracting the respective r values from the i values. The total 
areas and volumes are determined by combining the areas and volumes for 
each segment. For example 

v * = J? ( v ^ (ii • v «- (i) ) 

segment 

whe re 

V., = total trunk volume 

tk 


V tki (i * - 
V «kr<‘> = 
< V tkr<‘> = 


i value of trunk volume for ith segment 
r value of trunk volume for ith segment 
0 if ith segment is not in ground contact. ) 
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The trunk is symmetric about axes AA and BB (Fig. 2). 

Since roll motion is not considered and the ground profile is assumed to 
be two dimensional, only the right half of the trunk needs to be analyzed. 
The results of the left half will be similar. The right side of the trunk 
can be divided into four sections; two curved sections, each subtending 
90° at the center of curvature and two straight sections, on each side of 
axis AA. The curved sections are divided into N segments each, and the 
straight sections are divided into M segments each. Thus the complete 
trunk is divided into 4(N+M) segments. (See Figure 2. ) 

The following assumption is made in deriving the areas and 
volumes for each segment: 

The ground under any particular segment is con- 
sidered parallel to the hard surface, and at an 
elevation corresponding to the ground profile at 
the segment center projection (as shown in 
Figure 3). This assumption represents the ground 
surface and the hard structure of the cushion by a 
series of short, parallel sections which, when 
chosen sufficiently small, closely approximate the 
actual ground profile and hard surface orientation. 

In order to find the elevation of the segment center and the 
ground, the following five quantities are required. 

(a) The coordinates of the CG; Y , X 

eg eg 

(b) The pitch angle, $ 
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(c) The position of the CG with respect to the 
cushion center. 

(d) The distance of the segment center from the 

cushion center, X (i); i = 1, 2(M+N) 
cx 

(e) The ground profile, Y (i) as a function of X, (i). 

g IX 

C. 2. 2. 1 Hard Surface Clearance 

Figure 3 shows the various parameters involved 
in calculation of the hard surface clearance for each segment. 



For a given trunk orientation, X , 

Y cg and 0 

are known. 

From Figure 3, 



X cc = X cg 

+ GG Sin0 - CC Cos0 

(7) 


}H 

II 

}H 

- GG Cos0 - CC Sin^ 

(8) 


cc eg 


From Figure 2 



X (i) = - 

cx 

'L . 

-^-+ (| + R 2 Sin0 2 ) Cos6(i) 



for 1 < i < N 

(9) 

where 

fi(i) = (i-l)/3+ | 

(10) 

and 

= 90! 


(11) 
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( 12 ) 


where 


X (i) = 
cx ' 


- (i-l-N) del - del /2 

X Xi 


for N < i < M4N 


del = L /(2M) 
x s ' 

X (i) = (i-N-M-1) del + del /Z 

CX X X 

for M+N < i < N+2M 


and finally ^ 

X cx (i) = -j- + (| + R 2 Sin0 2 ) Sinfi (i-N-2M) 

for N+2M < i < 2N+2M 


and 6 (i) is given by Eq. (10). 


(13) 

(14) 


(15) 


Knowing X cc , Y cc and X cx (i), X h (i) and Y h (i) 
can be calculated as follows. 


x h (i) = 

X + X (i) Co 80 

cc cx ' 

(16) 

Y h (i) = 

Y cc + x «c<« Sia ° 

(17) 


i = 1, 2M+ 2N 



Since the ground profile is available in the 

following form, 


Y g = f(X g ) (18) 

substitution of the coordinates gives 

Y g(i) = f ( X h (i) + Y h {i) tan *) (19) 

i = 1,, 2M + 2N 
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ith segment is; 


Finally, the hard surface clearance for the 


Ygh^) = Y^ 1 )/ 0 ® 8 * - Y g (i) Cos 0 (20) 

i = 1, 2M+2N 

C. 2. 2. 2 Areas and Volumes 

C. 2. 2. 2. 1 Calculation of I Values 


As explained in Section C. 2. 2, 

i values of areas and volumes for a segment are calculated assuming that 
the segment is out of ground contact. The trunk sectors are shown in 
Figure 6a. From geometry, the trunk cross sectional area A tki<i> i* 
given by 


A tki^ A 1 " A 2 + A 3 " A 4 + A E 


where 


0 

A = ± R 2 
1 2 2 


A 2 = 


(R.-h ) 

R 2 Sin *2 


*,R' 


A 3 


1 1 


A 4 = 


Xb 


A 5 


(a-R 2 Sin0 z -X)(h y -R 1 ) 


and X = 


b(a-R 2 Sin0 2 ) 
b+h -6.. 

y i 


( 21 ) 
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position of the segment. 


The trunk volume depends on the 


and 




del 

x 

for 


A tki^ 

N < i < N+2M 


v adW = 




+ x e ) 


A tki<^ 


for i < N 
or i > N+2M 


where X is the horizontal distance of the centroid of the area 
e 

the inner trunk attachment point (X-coordinate of centroid). X 
lated as follows 


X = 
e 


A 1 X 1 ~ A 2 X 2 + A 3 X 3 ~ A 4 X 4 + A 5 X 5 
A tki (i > ~ 


where Xj, X^» etc. , are the X coordinates of the centroids of 

a a «i._ i.: i__ 

-Ti. J v • 9 1 


Xj = R 2 Sin 0 2 - 4Sin 2 (0 2 /Z) R £ /30 2 


X 2 = 0. 6667 R 2 Sin* 2 

X 3 = R 2 Sin0 2 + 4Sin 2 (0^2) R l /3* l 
X 4 = a - 0. 333X 

X = R Sin0 + 0. 333(a-R Sin0 -X) 


( 22 ) 


(23) 


A it . from 
tki 

is calcu- 


the areas 



given by 


The cushion area A c j^(i) is 


A chi (i) = <| + *- 2 S* V del x (24) 

for N < i < N+2M 

A chi (i » = <! + R 2 Si “»2> 2 I < 25 > 

for i < N 
or i > N+2M 


To calculate the trunk-to-cushion 

flow area, it is necessary to know the number of trunk holes inside the 
cushion. From Figures 1 and 6, the number of rows of holes communi- 
cating with the cushion is given by the integer value of (i^-fp/S^) + 1. 

The trunk-to-cushion flow area is thus given by 


A tkchi 


(i) 


integer 



+ 1 


N iA 


del 

x 

“3~~ 


(26) 


for N < i < N+2M 


where S = cushion periphery 

= 2L fl + 2ir(| + R 2 Sin0 2 ) 


(27) 


and 


A tkchi (l) 


= integer 


i_-i 

-ir 6 - + 1 


VS. 


0(^ + R 2 Sin0 2 ) 

s 


( 28 ) 
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for i < N 
or i > N+2M 


The trunk-to-atmosphere flow 


area is given by 


A tkati (i) = N r A k N h ' 


del 


" A tkchi (l) 


(29) 


for N < i < N+2M 


and 


A tkati (l) N r A h N h ' 


S(f +R 2 Sin* 2 ) 


A tkchi<‘> < 30 > 


for i < N 
or i > N+2M 


The cushion- to -atmosphere flow 


area is given by 


A gapi (l) = (Y gh (i)-h y ) del x 


(31) 


for N < i < N+2M 


where Y g ^(i) is calculated from Equation. (20), and 


A .(i) = (Y . (i)-h ) + R_Sin*,) 

gapi ’ 1 gh' 7 v 2 2 V 


gh 


(32) 


for i < N 
or i > N+2M 


given by 


Finally, the cushion volume is 


V chi (i) = Y gh (i) <7 + R 2 Sin V del x - (A r A 2 ) del x < 33 > 


for N < i < N+2M 
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and 


B , d 


v chi«> = Y g h (i) f <t + R 2 Si "V 


- /? ( 2 + X i2^ A l“ A 2^ 


( 34 ) 


where 
X 12 = 


X l A l” X 2 A 2 

_____ 


for i < N 
or i > N+2M. 


C. 2. 2. 2.2 Calculation of R Values 


The r values represent the changes 

in areas and volumes due to trunk-ground contact. Ground contact occurs 
when Y gh (i) < hy (see Figure 3). With ground contact, the r values are 
calculated as follows. The trunk cross sectional area i s given by 


A tkr^ A 6 “ A 7 + A 8 A( - 


(35) 


where A,, A_, etc. , are the areas of the sectors shown in Figure 6b. 
o 7 


A ~ 0 

a 6 2*3 


(WW’) 


R 2 Sin(» 3 


R^0 

l v 4 


. (vy-v 1 *) RiSin?4 
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and 


R - <h -Y .(i)) 

0, = cos -1 2 -jr g h — 

3 R 2 


R -(h -Y (i)) 
» 4 . Co,' 1 <■>- . 


The r value of the trunk volume 


V^ r (i) is calculated as follows 


W 1 * = A tkr^) del x 


( 36 ) 


for N < i < N+2M 


V tkr^ = *<! + X er> A tkr^ 


for i < N 
or i > N+2M 


where X is the X- coordinate of the centroid of area A., (i) 

er tkr' ' 


er 


A*X 4 - A„X„ +A 0 X 0 - A n X n 

w y l t v v / ? 

A., 

tkr 


and Xg, X ? , etc. , are X coordinates of the areas Ag, A^, etc. , respec- 
tively 


. 2 


X 6 = R 2 SinH> 2 - 4Sin (t^/2) R 2 /30g 


X ? = R 2 Sin?> 2 - 0. 333R 2 Sin0 3 


X g = R 2 Sin0 2 + 4Sin 2 (^/2) Rj/3^ 


X 9 = R 2 Sin0 2 + 0. 333RjSin0 4 
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given by 


The cushion area A c ^ r (i) is 


A chr (i) = del x R 2 Sin *3 


(37) 


for N < i < N+2M 


and 


A chr (l) = 2 ( ( 2 +R 2 Sin0 2 ) 


- (§ + R 2 Sin0 2 - R 2 Sin0 3 ) 


(38) 


for i < N 
or i > N+2M 

The number of trunk orifice rows 
communicating with the cushion is given by the integer value of 


f,-i - 0,R, 
2 p 3 2 


+ 1 


The r value of the area of the trunk orifices communicating with the 
cushion is given by 


A tkchr (i) = ' inte 8 er 


tkchi' 


jf.-i -0,R, 
-i- gP-1 * + 1 


N h A h 


del 


(39) 


for N < i < N+2M 


88 



and 


A tkchr<‘> ' 


tkchi 


(i) - integer 


W*3 R 2 


+ 1 


N lA 


+ R 2 Si “* 2 ) 

S 


( 40 ) 


for i < N 

i > N+2M 


Similarly, the number of orifice 

rows communicating with the atmosphere is given by the integer value of 


e r (*" ( 


V^r- 1 ^)) -+4 R 1 
£ 1 


+ 1 


The r value of the trunk- to- atmosphere area is thus 

V ('- S h)) 


A tkatr (i > = A tkati (i) ' int '* er 

del 


<l> R 
y 4^1 


♦"Vh 

(41) 


for N < i < N+2M 


and 


A tkatr {i) = A tkati (i > “ inte 8 er 




+ 4 R 1 


+ 1 


Vh 


+ R 2 Sin4> 2 ) 

s 


for i < N 
or i > N+2M 


( 42 ) 
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Equations (39) through (42) give 
the r values that correspond to a perfect seal during ground contact; i. e. , 
when the trunk orifices and clearance gap are completely blocked. How- 
ever, due to ground irregularities, trunk ribs and imperfections, brake 
pads, etc. , the seal will not be perfect and the actual r values will thus 
be smaller than those found from Equations (39) through (42). These 
reduced r values are related to the perfect seal values through blockage 
factors as follows. 

^tkchr^ ^tkchr^ 


x Per 


tk 


(43) 


perfect seal 


■^tkatr^ ^"tkatr^ 


x Per 


tk 


perfect seal 


(44) 


where Per^ is the blockage factor (less than unity) that depends on the 
trunk orifice sealing characteristic during ground contact. 


is given by 


The cushion clearance gap area 


given by 


A (i) = 
gapr' 


A .(i) 
gapi 


(45) 

The cushion volume r value is 


V chr <i » = - del * <A 6- A 7> ,46) 

for N < i < N+2M 

V chr<‘> ' ' ^<f + X cr" A 6- A 7» (47 > 
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for i < N 

i > N+2M 


where 


cr 


A 6 X 6 ‘ A 7 X 7 
A 6 - A 7 


area is given by 


Finally, the trunk-ground contact 


A tkcn (i) = ( R 2 Sin<<> 3 +R 1 Sin< * > 4 ) del x (48) 

for N < i < N+2M 


A tkcn* x) 2 (*2 +R 2 Sin<t> 2 + R 1 Sin 4 > 4 ) 

- (| +R 2 Sin4> 2 - R 2 Sin4> 3 ) 2 j (49) 

for i < N 
or i > N+2M 


The values of V^, A^, A tkch , 

A., V , and A for the full trunk and cushion are obtained by subtract- 
tkat ch gap 3 

ing the r values from the i values for each segment and summing them 
over all the segments. 


tk 


2(N+M) 

2 E 

i= 1 


V tki (i > 


V tkr< i 3 


(50) 


ch 


2(N+M) 

2 E 

i= 1 


A chi (i) 




(51) 


tkch 


2(N+M) 

2 E 

i= 1 


A tkchi^ A tkchr^ 


(52) 
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tkat 


2(N+M) 

* E 

i= 1 


A tkati (l) 


A tkatr (i) 


(53) 


V 


ch 


2(N+M) 

2 E 

i=l 


V chi< i ) - W 1 * 


+ V 


chd 


(54) 


where V ^ is the dead (inactive) cushion volume and 


2(N+M) 

A = 2 S 

gap i^l 


A .(i) - A (i) 
gapi' gapr 


(55) 


The factor 2 in the above equa- 
tions is included, because the expressions in brackets have been calculated 
for one-half of the (symmetrical) trunk. The r values are zero when 
ground contact does not occur, i. e. , Y gh (i) > h^. Otherwise they are 
calculated from Equations (35) through (47). 


Due to ground irregularities, trunk 

imperfections, brake pads, etc. , the cushion seal during ground contact will 

not be perfect (i. e. , A * 0). Therefore, a minimum value of A is 

gap g a P 

defined such that in ground contact, A = A^^. -^leak * S re ^ a ^ e< ^ to the 
equilibrium gap area through the cushion blockage factor, as shown below. 


A. . = Per A ' 
leak ch gape 


(56) 


where ^-g a p e the equilibrium gap area and Per^ i s the cushion blockage 


factor. 


C. 2. 2. 3 Center of Pressure 


The distances of the centers of pressure of 
each segment from the center of the cushion are required in order to 
estimate torques acting on the ACLS. The positions of the centers of 
pressure depend on whether the segment is in ground contact or not. 
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C. 2. 2. 3. 1 Segment Out of Ground Contact 


Such a situation exists for the 

ith segment when Y^(i) > h^. Since it is assumed that the pressures 
inside the trunk and the cushion are uniform, the pressure centers of a 
segment will coincide with respective centroids of the projected area. 


As can be seen from Figures 2 

and 4, the center of pressure for the ith segment of the cushion, X ch (i) is 

X ch (i> = - ~T “3? ( f +R 2 Sin4> 2 ) 

* Cos ( (i- 1 ) /3 + P/2) Sin(/3/2) (57) 

for i < N 


X ch< i > s X c*« 


for N < i < N+2M 


(58) 


and 

x ch (i) = + Tl (7 + R 2 S ^ 2 ) 

• Sin ( U-N-2M- 1)0 + J3/2) Sin(/3/2) (59) 


for i > N+2M 


C. 2. 2. 3. 2 Segment in Ground Contact 


Figure 4 shows the parameters 

needed to calculate the distances of the center of the cushion pressure X^U) 
and the center of trunk pressure X^i) for the ith segment in ground contact. 
From geometry. 
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for i < N, 



X ch (i) ' - -T - 3? ( l + V 1 ”*! - V**} 1 



Cos {(i-l)P + P/2) Sin( 0/2) 

( 60 ) 


L . 

X (i) = - - XX2 Cos ((i- 1)0 + 0/2) 

tit ^ 

(61) 

where 

_ 4 Sin(0/2) (RR 3 -RR1 3 ) 

a 2 2 

3 P (RR -RRl ) 



RR = j + R 2 Sin<t> 2 + R 1 Sin<t> 4 


and 

RRl = | +R 2 Sin<(> 2 - R 2 Sin<t> 3 



For N < i < N+2M 



X ch® = X c* (i) 

(62) 


X tk (i > = X cx (i > 

(63) 


Finally, for i > N+2M 



x ch (i) ' -T + Tfi <1 + R z Sin+ z - 



Sin ( (i-N-2M-l)/3 + 0/2) Sin(/3 /2) 

(64) 

and 

L, 

X (i) = + XX2 Sin((i-N-2M-l)/3 + 0/2) 

tit ^ 

(65) 

where XX2 is 

obtained from Eq. (61). 



i 

! 
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C.2. 3 Pressure -Flow- Force- Torque Relations 


The flow diagram of the cushion is shown in Figure C. 1. 
The flow through the upstream fan orifice is given by 

“fan = A atfn C af V ^ <«> 

where ^fan = v °l ume flow through the fan 

A a tfn = orifice area, atmosphere to fan inlet 
= fan inlet orifice discharge coefficient 
P a tf n “ fan entrance pressure (negative, gage) 

P = air density 

The fan pressure rise, P f , is given by 


P = P _ p 
fan plm atfn 


(67) 


where 


^plm = P^ enum P ress nre (gage). 


The fan flow is obtained from the fan characteristics, and 
is a function of the fan pressure rise. 


“fan = « P fan> < fc8 > 

The remaining flows are found as follows. 

“plat = A pla, C pa (69 > 
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where 


where 


where 


where 


where 


Q plat - bleed flow rate 

A plat ' bleed area 

C pa - bleed discharge coefficient 

Q fnpl " Q fan ' Q plat 

Q fnpl - fan -to -plenum flow rate 

Q fnpl ~ Q p ich + Q pltk 

Qplch ” Plenum-to-cushion flow rate 

Q pltk - plenum- to-trunk flow rate 


Q , , = A . , C 
plch plch pc 


i 


2( p -P ) 
plm ch 


P , - cushion ores sure laa.ae) 

CXI - — - • 


A plch " P lenum - t o-cushion orifice area 


n = a n \/ 2 ^ P plm“ P tk^ 
u pitk A P itk c P t V p — 

P tk - trunk pressure (gage) 

A pltk - plenum-to-trunk orifice area 


(74) 


where 


where 


determined 


where 


pltk 


tkch 


tkat 


^tkch + ^tkat 

trunk- to- cushion flow rate 

trunk- to- atmosphere flow rate 


^tkch "^tkch ^tc 


4 


2(P tk- P ch) 


(75) 


^tkat "^tkat ^ta 


4 


2 P 


tk 


(76) 


^chat ^plch + ^tkch 


(77) 


Q c hat “ cushion- to -atmosphere flow rate 


Q . = A C 

chat gap gap 




2 P 


ch 


(78) 


The static (vertical) force developed by the cushion is 
from a force balance. 


Forcn = <P ch A £h + P, k A^l Co.* 
P^ - cushion pressure, gage 


A , - cushion area = 

ch 

P fc k - trunk pressure. 


2(M+N) 

2 Xf ( A chi <i> 

1 = 1 ' 

gage 


(79) 


A chr (i) ) 


Afkcn - trunk-ground contact area 
<t> - pitch angle 


2(M+N) 

2 Xt A tkcn< i > 

1 = 1 
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From Figure 4, Torn, the torque developed by the cushion 
and trunk pressure, is given by 

2(M+N) 

Torn - E 2P ch (a^. 0) - A chr «) ) (x^fi) - cc) 


( A &c^)( 


+ 2P .k l A tkcn<‘Ul X tk (i > ' CC ) 


(80) 


where the i values of the areas, and the center of pressure distances are 
given by Equations (24), (25), (37), (38), (48), (49), and (57) through (65). 

Under equilibrium conditions, the total cushion force equals 
the aircraft weight, and the torque is given by the product of the weight and 
the distance between the CG and the geometric center of the cushion. Under 
this equilibrium loading, the aircraft orients itself at a particular X # 
and <|>. If ground profile and X are known, variables Y and <|> uniquely 
define the aircraft and ACLS position, and define the variables A j^U)* 

W 1 ’- Va 111 ' X ch«> and X tk« i »- Thus ’ 


and 


= f li < Y cg’ ♦ » 

(81) 

Achr' 1 * = f 2i < Y cg’ * » 

(82) 

Aiken' 1 ' = f 3i ' Y cg* ♦' 

(83) 

^h' 1 ' - f 4i < Y cg- +> 

(84) 

X .k«> = f 5i < Y cg’ 

(85) 


Thus, for the static solution, Equations (66) through (85) 
can be solved to determine the pressures, flows, areas, etc. , Qj an , 

Q fnpl’ Q plch’ Q pltk’ Q tkat’ Q tkch’ Q chat’ P atfn’ P fan’ P plm’ P tk' P ch’ 
Y cg’ A chi (i) ’ A chr (l) ’ A tkcn (l) ’ X c h (l) ’ X tk (l) ' 
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The heave stiffness can be found from the slope of the 
load-deflection characteristic. 


Stif 

ycg 


AForcn 



(86) 


where AForcn - change in normal force 

AY - change in CG elevation 
eg 

Similarly, the pitch stiffness is found from the torque - 
rotation characteristic. 


where 

and 


Stif 


phi 


ATorn 

ST 


(87) 


ATorn - change in torque due to change in location of load 
with respect to CG 
A<t> - change in pitch angle. 


C. 3 The Dynamic Model 

The dynamic behavior of the ACLS is determined from the simul- 
taneous solution of the equations describing the body dynamics and fluid 
mechanics of the cushion. 

C. 3. 1 Body Dynamics 

C. 3. 1. 1 Force Balance 


the ACLS consist of: 


During dynamic motion, the forces acting on 


(a) The cushion pressure force (P , A , ) Cos<f> 

ch ch 

(b) The trunk pressure force during trunk- 
ground contact (P^ -^ken^ < ^ os< l ) 

(c) The aircraft weight (M g) 

St 
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(d) The aerodynamic drag (1/2 • A , • p v 2 ) 

D ph * 9 

where is the heave drag coefficient, 

is the projected heave area and v is 

the heave velocity, dY /dt. 

eg 

(e) The trunk damping force during trunk- 
ground contact, Foret, given by 

2(M+N) 

Foret = 2 Xf Forct(i) 

i=l 

where 

Forct<i) = 4irafjjr + fir K* (i »- cc )) 

if segment is in ground contact. 

Forct(i) = 0, if segment is out of 
ground contact 


The basic equation of motion is then found from 
Newton's law as follows. 


dY 

M c 


* d, r 


^ = (p ch A ch + P tk A tkcn> Cos,|> - 


M a g 


- 1/2 C D A ph P 
C. 3. 1.2 Torque Balance 


Of- 


Foret 


(88) 


The torques acting about the CG of the aircraft consist of 
(a) The cushion pressure torque 

2(M+N) 


cp 


2 ^ P ch ( A chi (i) - A chr (i ») ( X ch (i »- CC ) 
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(b) 


The trunk pressure torque 


2 (M+N ) 

T t P = 2 (w ) ( x tk (i >- cc ) 


(c) 


The torque due to ground friction 


2(M+N) 


Torf = -2 


tk 


( A «kd (i > ) » ) 


where jx is the coefficient of friction between the trunk and the ground. 


(d) The torque due to aerodynamic drag force 

\2 


T df = 1/2 C D \h C 


Of 


'enf 


where C en ^ is the horizontal distance of the center of the aerodynamic drag 
force from the CG. 


where 


(e) 


The torque due to trunk damping, Torqt 


2(M+N) 

Torqt = 2 / . Torqt(i) 


= 4infNr (tF + ar ( x „ (i '- c c))( x t k ' i »- cc ) 


if the segment is in ground contact 


0, if the segment is not in ground contact 


A torque balance about the CG then gives 


j2 , 

Inert • — = T + T. + Torf + T - Torqt 
7 cp tp df ^ 


dt 


( 89 ) 


where Inert is the pitch moment of inertia about the CG. 
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C. 3. 2 Fluid Mechanics 


(a) Plenum 

The fluid system consists of three interconnected 
chambers; plenum, trunk and cushion; fluid resistances and a fan. From 
the polytropic pressure-density relation, 


(P . +P J 
plm at' 

4r~ 


= constant 


( 90 ) 


dP 




lm 


Taking time derivatives, 


K <P , + p a *) * P T 


plm at 7 K plm 
“plm 


dt 


( 91 ) 


Conservation of mass in the plenum requires that 


dt ( P plm V plm^ P fnpl ^fnpl " P plm ^plch " p plm ^pltk 

( 92 ) 

From Equations (91) and (92) 


d P . 

plm _ 

~dt 


K Pplm^at* 

P V 
plm plm 



( 93 ) 


Substituting p plm = P fnpl = p (mean density), the 
dynamic flow continuity equation for the plenum is as follows. 


d P 


Pi™ - 


dt 


K(P plm +P a.> 

V pta 


^fnpl “ ^plch 



( 94 ) 
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(b) 


Cushion 


The continuity equation for the cushion is 
similar to Equation (94), with an additional term to include the rate of 
change of cushion volume due to motion. 


d P 


It 


ch _ 


K(P ch +P at> 


^plch + ^tkch" ^chat 


d V 


ch 


at 


(95) 


d V 


where 


"eft" 


ch 


- rate of change of cushion volume. 


(c) 


Trunk 


dP tk. K < P tk +P at> 


dt 


tk 


d V 


Q pltk" Q tkch‘ Q tkat " ~3t 


tk 


(96) 


Equations (88), (89), (94), (95) and (96) define 
the dynamic heave model of the ACLS. The flows, areas, volumes and 
lengths needed to evaluate these equations are found from the relationships 
derived in Section C.2. 


C. 4 Analytical Simplifications 

The analysis described above has been developed for a general 
ACLS configuration. In using this analysis, several situations exist when 
the general relationsips can be simplified without loss of accuracy. These 
exist, for example, when some of the (user supplied) orifice areas are so 
large that the pressure drop across them is negligible, and need not be 
included in the computations. The program automatically determines the 
cases when such simplifications are possible, and modifies the basic 
analytical model accordingly to eliminate unnecessary computation. The 
modifications to the basic analysis are described below. 
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(a) When the cushion is high above the ground, the cushion 

pressure is equal to the atmospheric pressure. Although 
this result can be obtained from computation of Equation 
(78), significant computing reductions are achieved by 
assuming P ^ = 0 and dP^/dt = 0 (without computation) 
when the height is very large (i. e. , prior to the landing 
impact). Thus, at the beginning of the simulation, when 
the cushion- to- atmosphere gap is more than the ground 
effect gap*, the right-hand side of Equation (78) is not 
computed, but is set equal to zero. When the cushion 
enters the ground effect region (i, e. , gap less than 
ground effect gap), the above constraint on Equation (78) 
is removed. 


(b) 


When the trunk- to -plenum orifice is large - i. e. , the 
pressure drop is less than 2% of the upstream pressure - 
the plenum and trunk are treated as a single chamber, 
and Equations (94) and (96) are combined into a single 
equation. 


— P = — P = 
dt plm dt tk 


K (P plm tP 


<Vm +V tk> 


^tkch ^tkat 


at> / 


^fnpl " ^plch 


M 

IT-) 


(97) 


(c) During dynamic operation, when the sharp peaks in 

cushion pressure occur, the force on the trunk membrane 
at the cushion-trunk interface will reverse (i. e. , the com- 
putations of Equations (95) and (96) will indicate P ch > P tk‘ 
In such cases, trunk motion will tend to equalize the 
pressure difference, and this is included in the analysis 
by treating the trunk and cushion as a single chamber, 

^Gap at which the cushion pressure begins to increase above atmospheric 
pressure. This gap is specified by the user. 
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with P ch = P^, and combining Equations (95) and (96) 
into a single equation. 


Ip = 1 P 

dt tk dt ch 


K(P tk +P at> / 

V tk +V ch \ ^P ltk + Q P lch " °tkat 


_ Q _ y 

chat dt v ch 



(98) 


(d) For the duration of time when situations (b) and (c) above 

exist simultaneously, the plenum, trunk and cushion are 

treated as a single chamber (P , = P., = P , ) and 

' ch tk plm' 

Equations (94), (95) and (96) are combined into a single 
equation. 


I P =-I p a _d = _ 

dt plm dt tk dt ch (V 


K(P +P J 
plm at 


, +v i+V., ) 

plm ch tk 


( Q 


fnpl ^tkat “ ^chat 


1 V ch " 1 V tk <"> 



APPENDIX D - SUBROUTINE DESCRIPTIONS 


A flow chart of the overall program is shown in Figure A. 1 and 
discussed in Appendix A. Detailed descriptions of the main program and 
subroutines are given below. 

D. 1 Main Program 

The main program coordinates the static and dynamic simulation, 
carries out the data conversion, prints the static characteristics and deter- 
mines the initial conditions. Internally supplied data (discharge coefficients, 
etc. ) are read in directly. Other data are read in through subroutine 
PROGIO, which also prints out the data. 

The data conversion section of the program converts trunk orifice 
2 2 

area (AH) from in to ft , mass from lbs to slugs, initial pitch angle 
(PHU) from degrees to radians, initial angular velocity (DPHII) from 
degrees per second to radians per second, and atmospheric pressure (PAT) 
from psi to psf. Density (RHO) is calculated from the atmospheric tem- 
perature. IFAN is set equal to zero to start the program with unstalled 
fan operation. 

Subroutine TRUNK is called next to calculate the equilibrium shape 
of the trunk. ISHAPE is set equal to zero if the trunk configuration is 
infeasible and the program is terminated. Subroutine SEGMT divides the 
trunk into segments. Calculation of the distance of each segment center 
from the cushion center and assessment of the trunk and cushion areas 
and volumes (i values) associated with each segment is accomplished by 
calling subroutine SHAPE 1. 

The static characteristics of the ACLS are determined by calling 
subroutine FLOW1. ISAVE is set to zero if the configuration is con- 
sidered infeasible (i. e. , the static flow equations do not have feasible 
solutions). The equilibrium values of the variables (pressures, flows, 
etc. ) for the given ACLS loading, along with the values at different load 
levels and load positions are printed. 
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The variable IPP is set to zero if the difference between the 
equilibrium trunk and plenum pressure is less than 2%. In such a situation, 
the computations are carried out by modeling the plenum and trunk as a 
single chamber with uniform pressure. 

Initial conditions for the given landing configuration are either: 

(a) Supplied by the user, i. e. , XCGI, YCGI, PHII, VELXI, 
SINKRT and DPHI; or 

(b) Calculated by calling subroutines COORDN, PROFILE, 
CLRNCE and SHAPE2, e. g. , AGAP, VCH, VTK, etc. ; 
or 

(c) Estimated by interpretation from the static characteristics, 
i. e. , PCH, PTK and PPLM. 

The initial PCH is assumed to be 0 (psfg) if the initial gap area is more 
than the ground effect gap area. The parameter NQ, used in subroutine 
STEQU, controls the transition from the zero cushion pressure constraint 
as the ACLS moves into ground effect. 

Finally, the dynamic simulation. Subroutine DYSYS, is called to 
carry out the numerical integration for the state equations , and determine 
the pressures, flows, motion, etc. of the ACLS as a function of time. 

D. 2 SEGMNT 


Subroutine SEGMNT divides the trunk into 4(M+N) segments. The 
values of M and N are speciiied by the user. The straight section of the 
trunk is divided into 4M segments (see Figure 2), and the circular part of 
the trunk is divided into 4N segments. Since the trunk is symmetrical 
about axis BB, and since roll motion is not considered in this model, only 
the right half of the trunk is analyzed. The results for the left half of 
the trunk are found from a mirror image of the right half. 


108 


The distance of the center of each segment from the cushion 
center is calculated, using Equations (9) through (15), and stored in array 
XCX(I). 

D . 3 TRUNK 


Subroutine TRUNK calculates the trunk shape parameters Rl, R2, 
LI, L2, PHI1 and PHI2 from the user supplied parameters A, B, L and 
HY, as shown in Figure D. 1. 


Calculation of Rl, R2, LI, L2, PHI1 and PHI2 requires iteration. 
The iteration procedure is as follows: 


(i) 

(ii) 


Initial guess, R2 

-1 R2-HY 




+ HY 


PHI2 = Cos 


R2 


(iii) Rl = (A^R2_Sm(PH^- ± (g+HY]j 


/• , ™,„ _ ^ _-l Rl-B-HY 

(iv) PHI1 = Cos ' ky 


(v) L2 = L - PHI1*R1 

(vi) R2S = L2/PHI2 

R2S is compared with R2. If the difference is more than the 

R2+R2S 

tolerance RTOL, a new value of R2 (R2 = ^ ) is assumed and the 

procedure repeated until ABS(R2-R2S) < RTOL. For infeasible configura- 
tions, the iteration is stopped after 50 steps and the program is terminated. 
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D. 4 


SHAPE 1 


The flow diagram for the subroutine SHAPE 1 is shown in Figure 
D. 2. SHAPE1 calculates the peripheral length S and i-values of the areas 
and volumes (per segment) listed below. 

Trunk cross section area (ATKI(I)) 

Trunk volume (VTKI(I)) 

Cushion area (ACHI(I)) 

Trunk- to -cushion flow area (ATKCHI(I)) 
Trunk-to-atmosphere flow area (ATKATI(I)) . 

The above parameters are calculated from Equations (21) through 
(30) of Appendix C. The contact area ATKCNI(I) is set equal to zero 
since ground contact is not considered in SHAPE 1. 

D. 5 COORDN 


Subroutine COORDN calculates coordinates of the center of the 
cushion (XCC, YCC) from the coordinates of die CG (XCG, YCG) and the 
pitch angle (PHI) using Equations (7) and (8). The X and Y coordinates 
of the center of each segment are then calculated using Equations (16) 
and (17). 

D. 6 PROFILE 


Subroutine PROFILE contains the user supplied ground elevation 
coordinates. The subroutine then calculates the elevation of the ground 
at the projection of the segment center for each segment and stores it in 
the array YG(I). 

D. 7 CLRNCE 


Subroutine CLRNCE obtains the values of the Y coordinates of the 
segment centers, YH(I), from subroutine COORDN and the ground elevation 
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Figure D. 2. Flow Diagram of SHAPE 1 


3 




corresponding to each segment, YG(I), from subroutine PROFILE, and 
calculates the hard surface clearance for each segment, YGH(I), using 
Equation (20). 

D. 8 SHAPE 2 


Subroutine SHAPE2 determines various areas, volumes and center 
of pressure locations for each trunk segment from the hard surface clear- 
ance YGH(I). Those areas and volumes that have already been calculated 
in SHAPE 1 are not recalculated. 

The subroutine is divided into three parts, as shown in Figure D. 3. 

Part 1; In Part 1, i values of cushion volume VCHI(I) and cushion- 
to-atmosphere gap area AGAPI(l) are determined from Equations (31) through 
(34). These values have not been calculated in SHAPE 1, since they depend 
on trunk orientation, which varies with hard surface height. 

Part 2; In Part 2, initially it is determined whether ground- 
trunk contact for a particular segment has been made, i. e. , whether 
YGH(I) < HY. If ground contact has occurred, r-values of the areas and 
volumes are calculated, which account for the decrements in the area and 
volume parameters due to the contact. Equations (35) through (49) are 
used to calculate the r-values. If a segment is not in ground contact, the 
r-values are set equal to zero. In either case, the location of the centers 
of pressure (cushion and trunk) for a segment with respect to the cushion 
center, XCH(I) and XTK(I) respectively, are calculated from Equations (60) 
through (65). 

Part 3; In Part 3, the r-values of different areas and volumes 
of the segments are combined, and subtracted from the i-values, as shown 
in Equations (50) through (55). 
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DATA INPUT 



Figure D. 3. Flow Diagram of SHAPE2 
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Figure D. 3 (Continued). Flow Diagram of SHAPE2 





ATKCHF * 0. 0 


VTK. ACH, ATKH, 
ATKAT, VCH, AGAP, 
ATKCN *«t to 0. 0 
— — < -■- 
DO 30 1*1. 2(M>Nr 


VTK * VTK ♦ 2. 0*( VTKI(I) - VTKR(I)) 
ACH * ACH * 2 . 0*(ACHI(I) - ACHR(I)) 
ATKCHF - ATKCHl(I) - ATKCHR(I) 


ATKCHF < 0. i 


1 AKKCH * ATKCH ♦ 2. 0*ATKCHF 
ATKAT F = ATKAT 1(1) - ATKATR(I) 


^ ATKAT F < 0. Op 


ATKAT F * O.Ol 


(Contlnirtd next page) 


Figure D. 3 (Continued) Flow Diagram of SHAPE2 












ATKAT = ATKAT ♦ 2. O^ATKATF 
VCH * VCH 4 2. 0*(VCHI{I) - VCHR(I)) 

AGAP = AGAP ♦ 2. 0*(AGAPI(1J - AGAPR(I)) 
ATKCN * ATKCN 4 2. 0 - (ATKCN](I) - ATKCNR(I)) 



Figure D. 3 (Concluded 


Flow Diagram of SHAPE2 






FORCE 


D. 9 


The flow diagram of subroutine FORCE is shown in Figure D. 4. 

The following forces and torques are calculated in this subroutine. 

(a) FORCNs The normal pressure force is obtained by 
multiplying the cushion pressure and the trunk pressure 
by the respective areas and taking the vertical component. 

(b) TORN; The pressure force calculated for each segment 
is multiplied by the moment arm (i. e. , the distance 
between the centers of pressure and the CG). 

(c) FORCT: The r value of the contact area ATKCNR(I) 

is used to check whether ground contact for a 
particular segment has occurred. The trunk damping 
force is assumed to be equally divided amongst the 
trunk segments in contact. The damping force for each 
segment is calculated as shown in C. 3. 1. 1 (e). 

(d) TORQT: Torque generated by FORCT is calculated by 

multiplying the individual segment damping force by 
the center of pressure moment arm. 

(e) TORF: The contact friction torque is calculated by 

taking the product of the normal trunk pressure force, 
coefficient of trunk-ground friction, and vertical 
distance of the trunk contact zone from the CG. The 
variable TORQ is set equal to the sum of the torques 
TORN and TORF. 
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Figure D. 4. Flow Diagram of FORCE 











D. 10 FLOW1 


Subroutine FLOWT calculates the static performance characteristics 
of the ACLS. The height of the CG above the ground (YCG), pitch angle 
(PHI), cushion-to-atmosphere gap area (AGAP), and plenum, trunk and 
cushion pressures and flows are calculated for a range of loads and load 
CG positions including specifically the user supplied aircraft weight and 
CG position. 


D. 10. 1 Summary 

The subroutine essentially consists of three nested 
iteration loops, as shown in Figure D. 5. 

(i) Iteration 1; In Iteration 1, values of YCG 
and PHI are iterated until the ACLS load and 
torque equal (within a given tolerance) the 
aircraft weight and the static torque due to 
CG offset. 

(ii) Iteration 2; This iteration determines the 
various pressures and flows associated with 
the fan. plenum, trunk and cushion. It 
iterates the value of PFAN until the air gap 
required to satisfy the orifice pressure-flow 
relations is within a given tolerance of the 
air gap calculated by geometry. 

(iii) Iteration 3: This iteration is required (within 

Iteration 2) to determine the values of cushion 
pressure which satisfy the orifice pressure- 
flow relations. 
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Figure D. 5. Flow Diagram of FLOW1 


(Continued Mit page) 

















PSTART * PSTART ♦ PINC 1 
IR5T * IRST ♦ 1 
IRTT - IRST - 300 


| OCHAT - QPI.CH + QTKCHl 






Figure D. 5 (Continued). Flow Diagram of FLOW I 
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The subroutine begins by assuming a grid of YCG, PHI 
values to initialize Iteration 1. The values of pressures and contact areas 
are calculated for the assumed values of YCG and PHI using Iterations 2 
and 3. Subroutine FORCE is called next to determine the force and torque 
developed by the configuration tinder consideration. The CG offset corres- 
ponding to this torque is determined, and from this, a quadratic index of 
the difference between the ACLS force and aircraft weight, and the assumed 
and actual CG offset is formed. If the smallest quadratic index in the pre- 
vious grid is larger than the acceptable tolerance value, a new grid of 
YCG and PHI is formed and the procedure is repeated. Otherwise, the 
iteration is terminated after the calculations are repeated once more to 
obtain the final equilibrium values. The values of the performance vari- 
ables for an acceptable grid are stored along with the final equilibrium 
values and transferred to the Main Program for printout. 

D. 10. 2 Details 


D. 10. 2. 1 Iteration 1 


D. 10. 2. 1. 1 Initialization 


The initial YCG, PHI grid is 

formed as shown in Figure D. 6. The first value of YCG in this grid 
(YCGSTRT) is chosen equal to (HY+GG-0. 1). The last value (YCGSTOP) 
is chosen equal to (HY+GG+0. 1). Similarly, PHISTRT = -0. 1 (radians) 

and PHISTOP = 0. 1. An initial grid of nine points is thus formed, as 

shown in Figure D. 6. 


D. 10. 2. 1. 2 Calculations 


Subroutines COORDN, PROFILE, 

CLRNCE and SHAPE2 are called to determine the various areas and volumes 
associated with a particular grid point. Iterations 2 and 3 are then used 
to determine the various pressures and flows through the ACLS. If the 
configuration under consideration does not generate a positive value of 
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CG 



0 1st grid 
. % 2nd grid 

□ 3rd grid 
• Equilibrium value 

XDIS minimum for (3, 3) in 1st grid 

XDIS minimum for (2,2) in 2nd grid 

XDIS minimum for (2, 3) in 3rd grid and < XTOL. 

Figure D. 6. Grid Generation for FL.OW1 Iteration 
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plenum or cushion pressure for any feasible point on the fan curve, it 
is considered infeasible and the computations are carried out for the next 
grid point. For all feasible grid points, the normal force FORCN and 
torque about the center of the cushion, TORQ, are determined. The dis- 
tance of the CG from the center of the cushion, CC, is also calculated 
by dividing the torque by the force. FORCN and CC for a particular 
grid point are stored in two dimensional arrays, ZWT(IYCG, IPHI) and 
ZCC(IYCG, IPHI) respectively. For infeasible configurations, ZWT and 
ZCC are assigned a high value (= 10^) for subsequent detection and elimina- 
tion. The values of load, CC, YCG, PHI, AGAP, PCH, PTK and QFAN 
are also stored in arrays (static characteristics). 

D. 10. 2. 1. 3 Grid Point Evaluation 


After the arrays ZCC and ZWT 
are formed for all nine points of the grid, the points are tested to deter- 
mine which point comes closest to the actual values of ACES load and 
CG location (stored as CCI). 


A quadratic index is formed to 

determine the proximity of the calculated grid point from the physical value. 


XDIS = 


■ l 


\? 


f£WT(l, J) - L \ 

, ( ZCC(I, J) - CCI 1 

\ MASS *3 2. 2 / 

\ CCI / 


\ ? 


Calculation of XDLS is carried 

out for all the feasible grid points. From this, the minimum value of 
XDIS, corresponding to a particular grid point, designated (IS, JS) , is set 
equal to XTEST. If XTEST is within the iteration tolerance limit XTOL, 
the grid point (IS, JS) is taken as the equilibrium solution, and ICREST is 
assigned a value 2. The calculations are then repeated once more to 
determine updated values of YCG and PHI corresponding to the grid point 
(IS, JS), and the iteration is terminated. The values of YCG and PHI so 
found are returned to the Main Program. 
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If, however, XTEST > XTOL, 
a new grid is formed by observing the following rules. 

(a) If grid point (IS, JS) is not the midpoint, as in Grid 1 (Fig. D. 6), 
the grid is moved such that it does become a midpoint, as shown 
in Grid 2. 

(b) If the grid point (IS, JS) is the midpoint, the new grid is shrunk 
in size, so that each side is half that of the original grid and 
the midpoint remains unchanged. The transfer from Grid 2 to 
Grid 3 in Figure D. 6 illustrates this situation. 

(c) If none of the points in the original grid is feasible, a new grid 
is formed with higher values of YCGSTRT and YCGSTOP, but 
without changing the values of PHI. 

The procedure is then repeated 
as before until XD1S < XTOL. If after 40 iterations a solution cannot be 
found, the program terminates with an error message. 

D. 10. 2. 2 Iterations 2 and 3 


Iterations 2 and 3 find the value of PFAN such 
that the cushion-to-atmosphere gap area satisfying the orifice pressure- 
flow relations is equal (within a given tolerance) to the gap area generated 
by SHAPE2 in Iteration 1. 


The procedure is as follows. 

(a) Initially, two fan pressure rise increment values P1NC and PINC1 

are defined. PINC1 is used to converge rapidly to the first feas- 
ible configuration, starting from PFAN = 0. Thereafter, PINC is 
used to determine the static characteristics for 10 values of PFAN. 
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(b) 

(c) 

(d) 

(e) 

(f) 

This 

(g) 

(h) 

(i) 

(j) 

(k) 

( l ) 

(m) 

(n) 


QFAN is found from PFAN by calling subroutine CMPCRV (fan 
characteristic s ). 

PATFN is determined from Equation (66). 

PPLM is found from Equation (67). The value of PPLM is 
checked. If it is negative, the configuration is infeasible, 
PSTART is assumed to be PSTART + PINC1 and the process 
is restarted. 

QPLAT is found using Equation (69). 

QFNPL is found using Equation (70). 

Iteration 3 is set up next to evaluate PCH. 
is summarized in steps (g) through (n). 

Initially PCH is assumed to be equal to PPLM/2. 

QPLCH is found using Equation (72). 

QPLTK is found using Equation (71). 

PTK is found using Equation (73). 

QTKAT is found using Equation (76). 

QTKCH is found from Equation (74). 

PCHI is determined from Equation (75). 

PCHI found in (m) is compared with initially guessed PCH. If 

the difference is less than PTOL, the iteration is assumed to be 

PCH+PCHI 

complete. Otherwise, a new value of PCH = 2 i* 

assumed and the process repeated from step (h). 
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(o) If PCH is found to be negative, PST ART is assumed to be equal 
to PSTART 4- PINC1 and the process repeated from the beginning 
(Step (a)). 

(i) In case the new value of PSTART is just one step from 

PSTOP, a new value of PINC1 is defined as g-g 

so that the search for feasible configurations can be continued all 
the way to PSTOP. 

(ii) If a feasible configuration is not found within the full flow 
range of the fan, the calculations are performed for the next 
grid point. 

(p) If PCH > PTK, a new value of PSTOP is found, PSTOP = PFAN- 
PINC, and the procedure is repeated from the beginning (Step (a)). 

(q) QCHAT is found using Equation (77). 

(r) AGAPP(I) is found using Equation (78). 

(s) PFAN is stored as PFANS(I). 

The above procedure is repeated until 10 values 
of AGAPP(I) corresponding to 10 values of PFAN are found. Then the 
iteration value of AGAPP closest to the value of AGAP generated by SHAPE2 
is found. This value is designated as AGAPP(EDIF). If the difference 
between these values is less than ATOL, the iteration is complete; and 
updated values for the pressures, flows, and areas are calculated for the fan 
pressure value, PFAN(IDIF). If, however, | AGAPP(IDIF)-AGAP| > ATOL, 
new values of PFAN are chosen and the iteration is repeated from Step (a). 
The new initial value is set equal to PFANS(IDIF- 1), and the new final value 
is set equal to PFANS(IDIF+1). If, after 10 attempts, the iteration has not 
converged, the grid point in question is considered infeasible, and the cal- 
culations are then carried out for the next grid point. 


134 



D. 11 DYSYS 


Subroutine DYSYS, shown in Figure D. 7, coordinates the dynamic 
simulation and calls the various subroutines. Initially, DYSYS calls 
STEQU to set up initial values of the derivatives of the state variables. 
Then it calls RKDIF to get new values of the variables at the next time 
step. The values are printed after every MM (user supplied) time steps. 

DYSYS also carries out the following steps; 

(a) It determines the values of DVCH, DVTK, VELX and XCG. 

(b) It sets the value of the fan control parameter IF AN. 

This parameter is used by the fan subroutine CMPCRV 
to select the appropriate fan characteristic (IFAN = 0, 
unstalled operation; IFAN = 1, stalled operation). IFAN 
is set as follows 

QFAN < QP5 ; IFAN = 1 

QFAN > QP3 ; IFAN = 0 

(c) It sets the trunk inflation parameter IPCT, which deter- 
mines whether the cushion and trunk behave as two 
connected chambers or as a single independent chamber 
(See Appendix C). When PCH > PTK, the trunk membrane 
moves to equalize the trunk and cushion pressures. Thus 
for PCH > PTK, IPCT is set equal to zero and the compu- 
tation is carried out (in STEQU) by considering a single 
chamber for the (combined) trunk and cushion. When 
PCH < PTK, the trunk remains in its normal inflated 
shape. In this case, IPCT = 1. 
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D. 12 RKDIF 


RKD1F is the numerical integration subroutine which calculates 
the values of the state variables at time t+dt, given the values at time t, 
using a 4th order Runge Kutta method. The integration scheme is sum- 
marized below. 

(a) The iteration procedure starts with the values of the 
state variables y^, y^, etc. , at time t. 

y^t) i = 1, n 

(b) The slopes Dy^(t) are then determined from y.(t) by 
calling STEQU. 

E»y ± (t) = dy.(t)/dt 

dt 

(c) The values y. ^ at time t + — are then determined, 
y il = y i + Dy i * dt/2 

(d) The slopes Dy^j(t + dt/2) are then determined by 
calling STEQU and using the values of y^ found in 
(c) above. 

(e) The values at time t + dt/2 are then determined 
y i2 = y i + Dy il * dt/2 

(f) The slopes Dy.£ + dt/2) are then determined from 
STEQU using the values of y^ found in (e) above. 

(g) The values y.^ at time t+dt are then determined 

= n + D yi2 • dt 
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(h) The slopes Dy^ at time t + dt are then determined from 
STEQU using the values of found in (g) above. 

(i) Finally, the values of the state variables at time t + dt 
are found as follows 

Y i (t+dt) = y^t) + (Dy i + 2Dy u + 2Dy_ 2 + Dy i;} ) dt/6 

During each integration step (i. e. , to advance from t to t + dt), 
STEQU is needed four times to determine the slopes (b, d, f, g above). 
The fifth call for STEQU in DYSYS is to check whether PCH exceeds PTK. 

D. 13 STEQU 

The flow diagram for subroutine STEQU is shown in Figure D. 8. 
STEQU determines the values of the derivatives of the state variables by 
substituting the state variables in the state equations. Besides the seven 
state variables, PPLM, PCH, PTK, YCG, SINKRT, PHI and DPHI, the 
state equations need values of other variables (such as flows, forces, torques, 
areas and volumes) which are determined from the state variables. 

The areas and volumes for the given state variables are obtained 
by calling subroutines OOORDN, PROFILE. CLRNCE and SHAPE2. The 
flows are obtained from subroutine FL.OW2 and the forces and torques 
from subroutine FORCE. 

The appropriate dynamic equations are chosen depending on which 
one of the five conditions below prevails. 

(a) Normal operation. Equations (94), (95) and (96). 

(b) When the pressure drop across the trunk orifice is 
negligible « 2%), Equations (94) and (96) are replaced by 
Equation (97). 


139 




DATA INPUT 





4 





CALL COORDN 
CALL PROFILE 
CALL CLRNCE 
CALL SIIAPEZ 
CALL FLOW2 
CALL FORCE 




T 

/ ipp m Q \ 

F 


J 10 




4 ■■ 

DER(l) end DER(3) 
from (971 




DER(l) from (94) 
DER(3) from (96) 


QCHFT * QPLCII ♦ QTKCH - DVCH 
AD IF « (AGAP - GEC*AGAPE)«ID - 0.0001 



_ 1 ” 

IPP * 0 > 

4 

DER(I). DER(2) end DER(3) 
from (99) 


DER(2) end I)FR(3) 
from (98) 


T 

Figure D. 8. Flow Diagram of STEQU 


(Continued next pitge) 
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Figure D. 8 (Concluded 




Flow Diagram of STEQU 



(c) ACLS above ground effect zone - i. e. , PCH = 0, and 
DER(2) = dP ch /dt = 0. (Right hand side of Equation (95) 
set equal to zero. ) When the cushion enters the ground 
effect zone, the above constraints on cushion pressure 
are removed, and dP^/dt is changed in NQ steps from 
zero to the value found from Equation (95). 

(d) When PCH > PTK, the cushion and trunk are combined, 
and the value of IPCT is set to 0 in DYSYS. For 
IPCT = 0, Equations (95) and (96) are replaced by Equa- 
tion (98). 

(e) If both conditions (b) and (d) exist simultaneously, the 
plenum, trunk and cushion are combined and Equations 
(94), (95) and (96) are replaced by (99). 

Forces and torques due to aerodynamic drag are also calculated 
in STEQU and included in the respective state equations. 

D. 14 FLOW 2 

Subroutine FLOW2 calculates the flows through the various orifices 
from the pressures (PCH, PTK, PPLM) and the cushion-to-atmosphere 
gap area (AGAP), as shown in Figure D. 9. The subroutine solves the 
10 pressure-flow equations, Equations (66) through (70), (72), (73), (75), 
(76) and (78) for the 10 unknowns, QFAN, QPLAT, QFNPL, QPLCH, 
QPLTK, QTKAT, QTKCH, QCHAT, PATFN and PFAN. The solution pro- 
cedure is as follows; 

(a) QPLTK is calculated from Equation (73). 

(b) QPLCH is calculated from Equation (72). 

(c) QTKCH is calculated from Equation (75). 
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(d) QTKAT is calculated from Equation (76). 

(e) If the upstream orifice drop is not included in the 
model (i. e. , AATFN > 1 sq ft), PATFN = 0, 

PFAN = PPLM, and QFAN is determined from the 
fan subroutine CMPCPV. If the upstream resistance 
is included (AATFN < 1 ft^), an iteration is carried 
out to determine PATFN, as described in steps (f) 
through (i) below. 

(f) An initial value of PFAN is chosen and represented by 
the variable PEN(J). 

(g) PFAN is determined from Equation (67). 

(h) QFAN is determined from the fan subroutine CMPCRV. 

(i) A new value of PEN, PEN(J+1), is found using Equation 

(66). If the difference between PEN(J) and PEN(J+1) is 

less than PTOL, PATFN is set equal to PEN(J+1). 

Otherwise a new value, PEN(J+2) is determined and the 

procedure repeated from step (h). The new value 

PEN(J+ 2) is selected as ' ' 1 ^ ' ‘ • 11 tn ® sign 

of PEN(J)-PEN(J+1) is the same as the sign of 

PEN(J-2)-PEN(J-l). Otherwise PEN(J+2) is selected as 
PEN(J)+PEN(J-2) Thi , 

iteration scheme was chosen 
to bring about rapid convergence. 

(j) QPLAT is determined from Equation (69). 

(k) QFNPL is determined from Equation (70). 

(l) QCHAT is determined from Equation (78). 
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D. 15 CMPCRV 


Subroutine CMPCRV determines the fan flow (QFAN) for a given 
fan pressure rise (PFAN). 

The fan characteristics (see Figure 7) consist of two curves - 
the unstalled characteristic and a stalled characteristic. Both curves 
have been expressed in terms of 4th order polynomials. The polynomial 
coefficients and the pressures and flows indicating the curve limits of 
Figure 7 are supplied by the user. 

The region of operation is determined by QFAN. If QFAN is 
calculated to be less than QP5, IF AN is set equal to 1 and the fan 
operates in stall. If QFAN is more than QP3, IF AN is set to 0 and 
the fan operates in the unstalled region. If QFAN lies between QP3 
and QP5, IF AN is set equal to its previous value, and the fan continues 
to operate in the stall or unstalled region, depending on the region in 
which it was operating previously. IF AN is determined in DYSYS and 
is transferred to CMPCRV through COMMON. 
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ILLUSTRATIVE SIMULATION - 
INPUT DATA AND SAMPLE PRINTOUT 
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Landing Impact Simulation (Section 3, 2) 


Internal Data 

List 

CPA 

0. 6 

CAF 

1 . 0 

CPC 

0. 6 

CPT 

0. 9 

CTC 

0. 9 

CGAP 

1 . 0 

CTA 

0. 9 

CKK 

1.4 

GEC 

10. 0 

ZETA 

0. 05 

PERTK 

0. 85 

PERCH 

0. 15 
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0. 4 

DECCL 

15. 0 
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